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Abstract 


"\ 

The  navigation  performance  of  the  geometric  dilution 
of  precision  (GDOP)  GPS  satellite  selection  criterion  versus 
an  alternate  GPS  satellite  selection  criterion  for  an  inte¬ 
grated  GPS/INS  Navigator  in  a  jamming  environment  is  tested 
and  analyzed. 

The  theoretical  development  of  GDOP  and  GDOP-based  sat¬ 
ellite  selection  criterion  are  review  From  this  revie* 
a  satellite  selection  criterion  known  ^  the  cost  criterion 
(which  is  a  combination  of  noise-weig  c  GDOP  and  the  a 
priori  covariance  weighted  least  squares  position  error)  is 
selected  for  performance  testing  against  a  GDOP  satellite 
selection  criterion. 

Performance  assessment  of  the  two  satellite  selection 
criterion  is  accomplished  using  a  modified  version  of  the 
integrated  GPS/INS  (IGI)  Computer  Simulation.  The  IGI  Sim¬ 
ulator  is  modified  to  simulate  adaptive  band-width  control 
of  the  satellite  tracking  loops,  the  use  of  a  null  steering 
antenna,  the  use  of  an  alternate  satellite  selection  criter¬ 
ion,  and  the  effects  due  to  jammers.  The  performance  assess¬ 
ment  is  based  on  simulating  a  single  F-4  Close-Air-Support 
Mission  encountering  five  different  jamming  scenarios. 

v 

For  the  five  test  missions,  the  cost  criterion  demon¬ 
strates  a  performance  improvement  of  10-350  feet  in  mean 
radial  position  error  with  respect  to  GDOP  performance.  How¬ 
ever,  a  performance  degradation  of  10-30  feet  mean  radial 
position  error  is  observed  for  missions  with  rapidly  varying 
jamming  power. 


GLOBAL  POSITIONING  SYSTEM  SATELLITE 


SELECTION  EVALUATION  FOR  AIDED 
INERTIAL  NAVIGATION 

I  Introduction 

The  NAVSTAR  Global  Positioning  System  (GPS)  is  a  space- 
based  satellite  navigation  system  that  provides  highly  accu¬ 
rate  global  navigation  capability  to  suitably  equipped  users. 
The  GPS  User  Equipment  (UE)  set  collects  and  processes  GPS 
satellite  navigation  signals  to  determine  the  user's  three 
dimensional  position,  velocity  and  GPS  time.  The  GPS  UE  set 
can  operate  either  in  a  stand-alone  mode  or  an  integrated 
mode  (integrated  with  other  navigation  systems).  Integration 
of  the  GPS  UE  set  with  navigation  systems  such  as  Doppler 
radar  or  an  inertial  navigation  system  ^INS)  provides  accurate 
navigation  under  severe  operating  conditions  (i.e.,  high 
jamming  levels  with  respect  to  the  GPS  satellite  signals) 
(1:2).  The  benefits  of  an  integrated  GPS /INS  have  been  well 
documented  in  a  number  of  reports  (Refs  2,  3,  5,  6,  and  10). 

A  few  of  the  benefits  of  an  integrated  GPS /INS  navigator  are 
(3:145) : 

(1)  Increased  UE  set  tolerance  to  jamming 

(2)  Increased  UE  set  tolerance  to  severe  aircraft 
dynamics . 

(3)  Rapid  alignment  of  the  INS,  by  the  UE  set 


(4)  The  UE  set  updated  INS:  can  navigate  independently 
during  GPS  unavailability  and  provide  initialization  data  for 
rapid  GPS  reacquisition  when  GPS  signal  conditions  improve. 

The  concern  in  the  Air  Force  over  GPS  unavailability 
due  to  jamming  has  resulted  in  research  to  improve  the  GPS 
UE  sets'  anti-jam  capability  (Refs.  5,  6,  and  10).  This  re¬ 
search  resulted  in  a  number  of  methods  for  GPS  UE  set  improve 
ments.  From  the  number  of  methods  suggested,  three  main 
methods  for  improving  the  GPS  UE  sets'  jamming  tolerance  are 
(6:297) : 

(1)  adaptive  bandwidth  control,  inertial,  and/or  Doppler 
aiding  of  the  GPS  UE  sets '  satellite  tracking  loops 

(2)  adaptive  (pattern)  antennas  for  increasing  the 
strength  of  the  GPS  satellite  signal  and  decreasing  the 
strength  of  jamming  signals 

(3)  alternate  satellite  selection  criteria  to  choose 
satellites  which  will  maintain  high  navigation  accuracy  and 
not  be  jammed  off-the-air. 

All  three  methods  of  increasing  the  GPS  UE  sets'  capability 
in  a  jamming  environment  will  be  discussed  in  the  following 
pages.  The  primary  emphasis  of  this  report,  however,  will 
be  on  the  analysis  and  testing  of  an  alternate  satellite  se¬ 
lection  criterion  (to  be  known  as  the  "cost  criterion") 
against  the  normal  satellite  selection  criterion  also  used 
in  the  GPS  UE  set.  The  alternate  satellite  selection  crite¬ 
rion  versus  the  standard  satellite  selection  criterion,  will 
be  tested  using  the  rigorous  and  proven  integrated  GPS/INS 


(IGI)  computer  simulation  developed  at  the  Avionics  Labora¬ 


tory  at  the  Air  Force  Wright  Aeronautical  Laboratories, 
Wright-Patterson  AFB,  Ohio  (7). 

To  understand  the  GPS  UE  sets ’  satellite  selection  prob¬ 
lem,  a  thorough  review  of  the  terminology  and  theory  behind 
the  GPS  Program  must  be  accomplished.  Therefore,  the  GPS 
Program  will  be  reviewed  in  Chapter  II,  with  particular 
emphasis  on  the  composition  of  the  satellite  constellation, 
the  components  of  the  GPS  UE  set,  and  the  theory  and  opera¬ 
tion  of  the  UE  set.  The  satellite  selection  process  will  be 
further  isolated  and  analyzed  in  Chapter  III.  Chapter  III 
will  present  a  tutorial  on  the  classical  satellite  selection 
criterion  and  extensions  to  the  classical  criterion  to  improve 
integrated  GPS/INS  navigation  performance.  One  extension  to 
the  classical  satellite  selection  criterion  will  be  identified 
as  the  "cost  criterion",  which  will  be  tested  in  the  IGI  com¬ 
puter  simulation. 

An  overview  of  the  integrated  GPS/INS  (IGI)  computer 
simulation  will  be  presented  in  Chapter  IV.  This  chapter  will 
present  the  necessary  math  models  used  in  the  computer  simula¬ 
tion  as  well  as  other  information  necessary  to  follow  the 
data  analysis  in  Chapter  V. 

The  result  of  this  report  will  be  to  document  conclu¬ 
sive  results  of  the  performance  of  an  alternate  satellite 
selection  criterion  over  the  standard  satellite  selection 
criterion.  This  comparison  will  be  based  on  the  results  of 
a  rigorous  computer  simulation  using  a  realistic  environment 
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(a  tactical  aircraft  in  a  highly  dynamic,  jamming  environment) 
In  addition,  a  valuable  tool  will  have  been  created  (the  modi¬ 
fied  computer  simulation)  with  which  to  test  other  satellite 
selection  criteria  as  well  as  adaptive  receiver  and  adaptive 
antenna  designs. 


GPS  SYSTEM 


Figure  1.  GPS  System  Overview  (It 3) 


II  NAVSTAR  GPS  Program  Review 

Introduction 

The  NAVSTAR  Global  Positioning  System  (GPS)  is  a  space- 
based  satellite  navigation  system  that  provides  very  accurate 
three  dimensional  position,  velocity,  and  GPS  time,  world¬ 
wide,  to  suitably  equipped  users.  The  GPS  Program  is  com¬ 
posed  of  three  segments j  the  Space  Segment,  the  Control  Seg¬ 
ment,  and  the  User  Segment.  Each  segment  plays  an  important 
role  in  maintaining  the  24  hour-a-day,  all  weather,  high  ac¬ 
curacy  global  navigation  capability  that  is  NAVSTAR/GPS 
(Figure  1).  The  following  sections  in  this  chapter  will  dis¬ 
cuss  each  segment,  with  primary  emphasis  on  the  Satellite 
and  User  Segments.  The  level  of  detail  that  is  covered  in 
each  section  is  determined  by  the  further  development  of  ter¬ 
minology  and  theory  in  the  next  chapter. 

Space  Segment 

The  Space  Segment  is  responsible  for  planning  and  de¬ 
ployment  of  the  18  satellites  in  the  GPS  constellation.  The 
18  GPS  satellites  occupy  six  different  orbital  planes  with 
three  satellites  evenly  spaced  within  each  orbital  plane. 

The  six  orbital  planes  are  inclined  at  an  angle  of  55°  and 
each  is  separated  by  60°  of  longitude.  The  relative  phasing 
of  the  satellites  from  one  orbit  plane  to  the  next  is  40°  in 
an  easterly  direction.  The  phasing  of  satellites  means  that 
when  a  satellite  in  one  orbital  plane  crosses  the  equator, 
the  satellite  in  the  next  orbital  plane  (moving  in  an 


easterly  direction)  is  North  of  the  equator  by  40°  latitude. 
Each  satellite  operates  at  an  altitude  of  10,900  nautical 
miles  in  a  circular  orbit  with  a  period  of  one  satellite 
revolution  around  the  earth  equal  to  approximately  12  hours 
(Figure  2) .  This  particular  arrangement  of  satellites  in  the 
GPS  constellation  was  chosen  for  overall  user  accuracy,  the 
minimum  size  and  duration  of  regions  of  poor  navigation  per¬ 
formance,  and  survivability  of  the  satellite  constellation. 
Other  factors  favorable  to  this  constellation  are  the  ease 
of  replenishment  and  sparing,  ease  of  build-up  to  the  full 
constellation,  and  the  constellation’s  potential  to  grow  to 
24  satellites  (should  funding  permit)  (Ref  4  and  8). 

The  GPS  satellite  navigation  signal  uses  the  communica¬ 
tion  techniques  of  spread  spectrum,  pseudo-random  noise  (PRN) 
and  biphase-shift -Keying  (BPSK) .  A  spread  spectrum  signal 
will  typically  occupy  a  signal  bandwidth  much  larger  than 
the  information  bandwidth  (typically  greater  than  10  times 
the  information  bandwidth) .  Spreading  the  signal  has  no  ad¬ 
verse  effect  on  performance  but  provides  benefits  in  signal 
security,  anti- jamming  capability  (against  impulse,  contin¬ 
uous,  burst,  swept,  narrowband,  and  wideband  noise),  and 
transmission  security  (i.e.,  beyond  some  range  from  the  sat¬ 
ellite  the  transmitter  signal  is  buried  in  natural  back¬ 
ground  noise).  Pseudo-random-noise  (PRN)  and  biphase-shift¬ 
keying  (BPSK)  of  the  carrier  signal  are  analogous  to  time 
tagging  the  satellite  signal.  The  phase  of  the  carrier  sig¬ 
nal  is  periodically  shifted  forward  or  backward  as  determined 


by  the  instantaneous  value  of  a  very  long  sequence  of  digital 
ones  and  zeros.  This  sequence  is  called  PRN  because,  to  the 
casual  observer,  the  ones  and  zeros  appear  to  occur  in  a 
purely  random  fashion.  In  actuality,  the  code  generated  is 
predictable  relative  to  the  time  it  was  started. 

The  user  can  deduce  when  the  code  was  transmitted  by 
matching  his  own  copy  of  the  code  to  the  incoming  signal. 

The  amount  the  user  must  shift  the  code  to  match  the  incoming 
signal  determines  the  estimate  of  the  time  that  the  signal 
took  to  reach  the  user.  A  range  measurement  (from  the  user 
to  the  satellite)  can  be  calculated  by  multiplying  the  esti¬ 
mate  of  the  satellite  signal  propagation  time  to  the  velocity 
of  an  electromagnetic  wave  (approximately  equal  to  the  speed 
of  light) . 

The  GPS  satellites  transmit  navigation  signals  on  two 
frequencies,  1575.42  MHZ(Ll)  and  1227.6  MHZ(L2).  The  reason 
for  transmission  on  two  frequencies  is  that  during  daylight 
hours  solar  radiation  produces  a  belt  of  ionized  particles 
in  a  portion  of  the  atmosphere  known  as  the  ionosphere  (40- 
300  miles  above  the  earth  surface).  Signals  passing  through 
this  region  are  refracted  resulting  in  longer  than  normal 
propagation  of  the  satellite  signals  to  the  user  (and  a  re¬ 
sultant  longer  time  delay  estimate) .  The  increase  in  time 
delays  translate  into  ranging  errors  which,  if  uncorrected, 
can  sometimes  lead  to  relatively  large  position  errors.  The 
ionospheric  effect  has  a  predictable,  annual  pattern  of  var¬ 
iation.  However,  the  effect  is  not  a  totally  predictable 


phenomena.  For  high  accuracy  positioning  the  phenomena  cannot 
be  completely  modeled.  By  making  range  measurements  on  both 
LI  and  L2  signals  at  the  same  time,  a  simple  mathematical 
technique  can  be  used  to  correct  the  ionospheric  error 
(22i 115) . 


At  LI  : 


Ri  =  rtrue  +  (k/f  1  * 


At  L2» 


R2  =  RTRUE  +  ^k//f2  ^ 


Combining  and  rearranging: 


TRUE 


_  R1-(f2/f1)‘  R2 

l-(f2/fl)2 


Where:  and  R2  are  estimates  of  the  true  range,  R^^ 

fj  and  f2  are  the  respective  LI  and  L2  frequencies 
k  is  a  constant  associated  with  the  ionospheric  error 


Superimposed  on  the  LI  frequency  are  two  uniquely  coded 
signals,  a  precise  or  P  code,  and  a  coarse/acquisition  or  C/A 
code.  The  P  code  is  a  very  long  sequence  of  digital  pulses 
which  does  not  repeat  for  approximately  280  days  (a  PR N  code). 
The  pulse  train  is  created  by  a  complex  set  of  equipment  on 
the  satellite.  The  code  is  generated  at  10.23  megabits  per 
second  chipping  rate,  where  a  chip  is  the  time  interval  of 
one  pulse  in  the  pulse  train  (a  chip  for  the  P  code  equals 
97.7  x  10  seconds).  A  unique  portion  of  this  code  (one 
week)  is  assigned  to  each  of  the  GPS  satellites.  Therefore, 
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the  user  can  distinguish  from  the  code  which  of  the  satellites 
it  is  tracking.  The  C/A  code  is  a  short  sequence  relative 
to  the  P  code  and  has  a  chipping  rate  of  one  tenth  the  P  code. 
The  C/A  codes  are  chosen  from  a  family  of  distinct  codes 
called  gold  codes.  The  use  of  gold  codes  assures  minimum 
interference  between  satellites,  again  allowing  unique  satel¬ 
lite  identification  (Ref.  23).  The  C/A  code  was  chosen  to 
assist  the  GPS  UE  set  in  reducing  the  time  necessary  to  ac¬ 
quire  and  track  the  longer  P  code  (Figure  3)  (Ref.  23). 

Since  the  P  and  C/A  codes  are  chosen  not  to  interfere 
with  each  other  (i.e.,  minimum  correlation  among  the  codes), 
they  can  be  modulated  onto  the  same  carrier  frequencies  in 
the  satellite.  Transmitters  are  known  to  be  much  more  ef¬ 
ficient  with  constant  amplitude  signals.  Therefore,  the  P 
and  C/A  carriers,  though  derived  from  the  same  source,  are 
phase  shifted  90°  apart,  modulated  by  the  P  and  C/A  codes 
respectively  and  then  combined.  This  process  is  known  as 
phased  quadrature,  and  it  produces  a  composite  continuous 
wave  (CW)  signal  at  1575.42  MHZ  (LI)  (See  Figure  4). 

Both  the  P  and  C/A  codes  are  modulated  by  a  50  bit-per- 
second  data  message  by  the  method  of  biphase-shift-keying 
(BPSK).  The  total  block  of  data  transmitted  is  1500  bits, 
thus  the  user  will  require  at  least  30  seconds  to  read  the 
full  block  of  data  for  each  satellite.  The  data  message  con¬ 
tains  system  time,  satellite  clock  offset  and  drift  errors, 
satellite  position  information  (ephemerides ) ,  C/A  to  P  code 
handover  information,  and  an  almanac  of  the  health  and 


C/A  CODE  . . .  EASY  TO  ACQUIRE  BECAUSE 

LENGTH  *  1023  CHIPS 


THESE  STATES  INHERENTLY  SYNCHRONIZED  BY 
GENERATION  PROCESS  ABOARD  THE 
SATELLITE.  THIS  FACT  IS  KNOWN 


TO  GPS  USER. 


P  CODE  . . .  DIFFICULT  TO  ACQUIRE  BECAUSE 

LENGTH  ESSENTIALLY  INFINITE 


Figure  3.  Relationship  Between  GPS  Coarse/Acquisition 
Code  and  the  GPS  Precision  Code  (23:6). 


Signal  Formation  (23:10) 


positions  of  all  other  satellites  in  the  GPS  constellation. 
Relative  to  the  ranging  codes  (P  and  C/A  codes),  the  data 
message  rate  is  very  slow.  It  can  thus  be  superimposed  on 
the  ranging  codes  (and  separated  at  the  user’s  receiver) 
without  affecting  either  the  ranging  codes  or  the  data  mes¬ 
sage.  The  particular  ranging  code  and  data  message  are  com¬ 
bined  prior  to  modulating  the  carrier.  The  combination  of 
the  two  signals  is  by  modulo-2  addition  (exclusive -OR) .  For 
convenience,  the  same  set  of  data  that  is  modulo-2  added  to 
the  P  code  is  also  added  to  the  C/A  code.  The  generation  of 
the  LI  signal  (1575.42  MHZ)  containing  both  P  code  and  the 
C/A  code  with  the  data  message  added  to  both  has  just  been 
described.  The  L2  signal  (1227.6  MHZ)  is  created  in  the 
same  manner,  but  contains  the  P  code  and  data  message  only 
(no  C/A  code) . 

Control  Segment 

The  Control  Segment  is  composed  of  seven  Monitor  Sta¬ 
tions  (MS),  three  Ground  Antennas  (GA) ,  and  a  Master  Control 
Station  (MCS).  The  Monitor  Stations  and  Ground  Antennas  are 
stationed  throughout  the  world  on  U.S.  controlled  property. 
The  Master  Control  Station  is  located  at  the  Consolidated 
Space  Operations  Center  (CSOC),  outside  Colorado  Springs, 
Colorado. 

The  Monitor  Stations  continuously  collect  navigation 


signals  from  all  the  satellites.  This  information  is  then 
transmitted  back  to  the  MCS  for  processing.  At  the  MCS,  the 
satellites'  corrected  positions  and  clock  error/offset 


information  are  determined.  This  information  is  relayed  to 
a  Ground  Antenna  for  transmission  to  the  satellites. 


In  addition  to  the  corrected  navigation  data  the  MCS 
also  performs  satellite  station  keeping,  health  status  checks 
of  the  satellite  hardware,  and  many  other  satellite  mainten¬ 
ance  functions  required  to  keep  the  satellites  working  ef¬ 
ficiently.  This  satellite  information  is  transmitted  and 
received  on  a  coded  S-band  frequency. 

User  Segment 

The  purpose  of  the  User  Segment  is  to  develop,  test, 
and  produce  the  User  Equipment  (UE)  set.  The  UE  set  is  com¬ 
posed  of  the  control/display  unit  (CDU),  the  receiver/proces¬ 
sor  unit  (RPU) ,  and  either  a  fixed  reception  pattern  antenna 
(FRPA)  or  a  controlled  reception  pattern  antenna  (CRPA) ,  with 
an  antenna  control  unit. 

User  Equipment  (UE)  Operation:  Theory .  The  UE  set  can 
locate  its  three  dimensional  position  by  knowing  the  distance 
or  range  from  three  reference  points.  In  GPS,  the  reference 
points  at  any  time  are  positions  of  GPS  satellites.  A  UE 
set  can  deduce  the  range  to  a  satellite  at  any  time,  because, 
as  stated  in  the  satellite  signal  description,  the  time  of 
transmission  of  the  signal  can  be  determined  by  the  code. 

By  noting  the  time  of  arrival  of  the  signal,  the  user  can 
calculate  the  range  (See  Figure  5). 
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WHICH  IS  TRANSMITTED  . . . 


DEMODULATED . . . 


PRODUCING  THE  SAME  CODE 
AT  THE  USER.  BUT  DELAYED  . . . 


TIME  RECEIVED 
(AS  OBSERVED) 
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TIME  TRANSMITTED 
IAS  DEDUCED  FROM  CODE  STATE) 

Figure  5.  GPS  User  Equipment  Operation  (23:14) 
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Range  (R)  =  Speed  of  Light  (c)  x  (time  of  arrival 

-  time  transmitted)  (4) 

Equation  (1),  however,  assumes  that  the  satellite's 
clock  and  the  UE  clock  are  synchronized.  Therefore,  the 
calculated  range  will  be  in  error  by  an  amount  proportional 
to  the  time  error  (or  time  bias)  due  to  the  mismatch  of  the 
satellite  and  receiver  clocks.  Thus, 

PR  =  RT  +  RB  =  c  x  (ta-tT+tB)  (5) 

RB  =  c  x  tB  =  time  bias  ranging  error  (6) 

Where 

PR  =  Pseudo-range  measurements 
Rt  =  True  Range 

R_  =  Range  Error  due  to  time  bias 
c  =  Speed  of  light  in  a  vacuum 
t  =  Time  of  signal  arrival 

a 

t,p  =  Time  of  signal  transmission 
tg  =  Time  bias 

The  time  bias  introduces  another  unknown  into  the  solution 
of  the  range  equations,  in  addition  to  the  three  desired  com¬ 
ponents  of  the  user's  position.  To  allow  instantaneous  cal¬ 
culation  of  these  four  navigation  unknowns,  four  independent 
measurements  of  range  are  needed  to  four  satellites  (Figure  6) 
Because  of  the  time  bias,  the  range  measurements  are  known 
as  pseudoranges. 


WITH  AN  UNSYNCHRONIZED  CLOCK.  THE  GPS  USER  MUST 
TAKE  PSEUOORANGE  MEASUREMENTS  TO  FOUR  SATELLITES. 

*1*R1  +  Ctb 

■  R2  ♦ 

R3  ■  R3  ♦  ««b 
R4  •  R4  ♦  ctb 


Figure  6.  GPS  User  Equipment  Operation* 

Pseudorange  Measurements  (1*15) 


The  receiver  will  collect  and  calculate  pseudorange 
measurements  (simultaneously  or  sequentially)  to  four  selected 
satellites  as  well  as  changes  to  the  pseudoranges  as  a  func¬ 
tion  of  time.  This  latter  measurement  is  called  a  delta 
pseudorange  (or  simply  a  delta  range)  and  may  be  calculated 
by  counting  the  number  of  carrier  cycles  that  occur  during  a 
finite  interval  (usually  0.1  to  0.5  seconds)  or  by  differen¬ 
cing  pseudorange  measurements  at  the  beginning  and  end  of  a 
finite  time  interval. 

To  make  pseudorange  measurements,  the  receiver  must 
lock-on  and  track  the  incoming  signals,  determine  the  time-of- 
arrival  and  time-transmitted  parameters  for  pseudorange  cal¬ 
culation.  To  accomplish  these  tasks,  the  UE  essentially 
reverses  the  modulation  process  of  the  satellite  transmission. 
The  spread  spectrum  is  collapsed  back  to  a  single  frequency, 
continuous  wave  signal.  The  receiver/processor  unit  then  uses 
correlation  techniques  to  track  the  satellite  code  and  car¬ 
rier.  A  phase  lock  loop  is  used  to  track  the  carrier  fre¬ 
quency,  while  a  code  tracking  loop  (an  internal)  code  gener¬ 
ated  to  match  the  transmitted  satellite  code)  is  used  to 
track  the  coded  signal  (See  Figure  7).  The  receiver  detects 
maximum  correlation  between  the  two  codes  when  the  internal 
code  generator  is  perfectly  aligned  with  the  received  code. 
Misalignment  of  the  codes  causes  the  correlation  of  the  two 
codes  to  decrease  to  a  relatively  small  value.  The  preci¬ 
sion  with  which  the  time-of -arrival  can  be  measured  depends 
on  the  chip  size  (pulse  width)  of  the  code.  The  P  code  has 


RANGE  DETERMINATION 


Figure  7.  UE  Set  Operation i  Range  Deterknination/teeceiver  Operation  (1*31) 


a  chip  size  of  ±  97.7  nanoseconds,  which  is  equivalent  to 
±  29.3  meters  ranging  error.  As  stated  earlier,  synchro¬ 
nizing  on  the  P  code  is  very  difficult  due  to  the  fact  that 
the  P  code  takes  so  long  to  repeat.  The  C/A  code  repeats 
every  millisecond  and  has  a  pulse  width  ten  times  the  P  code, 
thus  its  primary  use  is  for  fast  acquisition  of  the  satel¬ 
lite  signals. 

In  general,  prior  to  tracking,  the  signal  generated  in 
the  code  tracking  loop  will  not  correlate  with  the  received 
satellite  signal  due  to  the  unknown  time  delay  for  the  satel¬ 
lite  signal  to  reach  the  user  and  the  unknown  UE  sets'  clock 
bias  errors.  Similarly,  the  signal  generated  in  the  phase 
lock  loop  differs  from  the  received  signal  due  to  an  unknown 
Doppler  offset  of  the  carrier  frequency  caused  by  the  rela¬ 
tive  velocity  between  the  user  and  the  satellite  and  a  bias 
in  the  UE  sets'  frequency  standard  (li30).  During  satellite 
acquisition,  the  UE  set  shifts  a  code  and  carrier  search 
window  (the  size  of  the  search  window  is  determined  by  how 
well  the  satellites'  positions  and  velocities,  and  the  user's 
position  and  velocity  are  known)  until  a  received  signal 
matches  a  known  satellite  signal.  Matching  a  satellite  sig¬ 
nal  in  the  acquisition  phase  results  in  an  estimate  of  the 
time  and  frequency  bias  errors,  time  delay,  and  frequency 
shift  of  the  received  satellite  signal.  Satellite  tracking 
uses  the  delays,  bias  errors,  and  offsets  as  measurement 
variables  from  which  user  position,  velocity,  and  GPS  time 


are  determined.  The  UE  set  accomplishes  satellite  tracking 
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using  two  distinct  modes  of  operation:  coherent  and  non¬ 
coherent  tracking.  The  following  description  of  satellite 
tracking  is  just  one  method  of  implementation. 

The  coherent  mode  of  satellite  tracking  uses  a  tau- 
dithered  tracking  loop  where  the  code  phase  is  alternately 
advanced  and  retarded  by  half  a  code  chip.  The  sharpness  of 
the  correlation  pulse  (Figure  8)  and  the  size  of  the  code 
chip  determines  how  accurately  the  received  signal  can  be 
matched  to  the  UE  sets'  copy  of  the  satellite  signal.  The 
dithering  of  the  correlation  pulse  permits  the  derivation  of 
an  error  signal  to  be  used  by  the  code  loop  to  maximize 
signal  correlation  (Figure  8) .  The  coherent  mode  uses  one 
correlator  to  do  tau-dithering  and  uses  the  second  correlator 
to  remove  the  code  from  the  signal  for  carrier  tracking  (each 
satellite  tracking  loop  has  two  correlators)  (5:15-18). 

The  code  clock  and  carrier  frequencies  of  the  satel¬ 
lite  signals  are  precisely  related,  L2  is  exactly  120  times 
the  clock  frequency.  Thus,  when  the  phase  lock  loop  is  track 
ing  the  carrier,  the  code  clock  is  further  inherently  phase 
locked  to  the  received  code.  This  configuration  allows  the 
carrier  loop  to  track  the  user’s  dynamics,  reducing  the  dy¬ 
namics  the  code  loop  must  track. 

The  coherent  mode  is  the  normal,  preferred  mode  of 
operation.  An  alternate  mode  of  operation  occurs  when  the 
carrier  loop  breaks  lock  due  to  severe  user  dynamics  and/or 
low  signal-to-noise  (S/N)  ratio  (primarily  due  to  jamming 
noise).  In  this  alternate  mode,  called  the  non-coherent 
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Sync  tracking  by  early-late  correlators 


Figure  8.  Tau-Dither  Tracking 
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mode,  the  UE  set  uses  a  second  order  code  tracking  loop. 

The  second  order  code  loop  uses  simultaneous  early-late  cor¬ 
relators  to  track  the  code.  The  second  order  code  loop  pro¬ 
vides  much  more  accurate  pseudorange  measurements,  but  must 
also  track-out  all  residual  signal  dynamics  and  jamming  noise 
to  maintain  code  lock.  The  bandwidth  of  each  code  loop 
greatly  affects  the  performance  of  the  loops.  The  bandwidth 
design  requires  compromise  between  the  conflict  of  designing 
to  track  maximum  dynamics  (large  bandwidth)  or  to  overcome 
increased  noise  due  to  jamming  (small  bandwidth) .  Early  UE 
set  designs  used  fixed  bandwidths  based  on  the  dynamics  and 
jamming  levels  expected  for  particular  groups  of  users  (low, 
medium,  or  high  dynamics  with  low  to  moderate  jamming) .  More 
recent  designs  have  employed  adaptive  bandwidth  features  (Refs 
5,  6,  and  10).  The  specific  bandwidths  in  an  adaptive  re¬ 
ceiver  are  determined  by  the  estimate  (or  measurement)  of 

the  users  dynamics  and/or  the  jamming  noise.  In  this  way 

« 

the  UE  set  can  optimize  its  operation  for  the  changing  envi¬ 
ronment  and  increase  the  set’s  navigation  accuracy  (5:14-21). 

When  a  GPS  set  is  integrated  with  an  aircraft’s  INS,  a 
number  of  positive  changes  result  (as  already  mentioned  in 
Chapter  I).  The  INS  information  provides  an  accurate  esti¬ 
mate  of  position  and  velocity  for  prepositioning  of  the  car¬ 
rier  and  code  loops  for  acquisition  of  the  satellites.  The 
code  loop  uses  the  position  estimate  to  determine  an  approx¬ 
imate  range  to  the  satellite  (or  more  specifically,  a  time 
delay  between  an  incoming  satellite  code  and  the  user’s 
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version  of  the  code)  from  which  to  start  searching  for  the 
incoming  signal.  The  carrier  loop  uses  the  velocity  estimate 
to  calculate  an  approximate  Doppler  frequency  shift  from 
which  to  start  searching  for  the  incoming  carrier  frequency. 
The  INS  aiding  of  the  code  and  carrier  loops  as  described 
above  result  in  faster  acquisition  and  tracking  of  the  satel¬ 
lites.  Once  tracking,  the  INS  continues  to  provide  informa¬ 
tion  on  the  aircraft’s  dynamics  to  the  carrier  loop,  thereby 
reducing  the  amount  of  dynamics  the  carrier  loop  will  have 
to  track.  Therefore,  the  UE  sets*  carrier  loop  can  use 
smaller  bandwidths  and  still  maintain  lock  during  severe  air¬ 
craft  dynamics.  In  the  coherent  mode  (when  the  coupled  car¬ 
rier  loop-INS  is  tracking),  the  code  loop  will  benefit  from 
the  increased  capability  of  carrier  tracking.  In  the  non¬ 
coherent  mode  (when  the  carrier  loop  can  no  longer  maintain 
track  even  with  INS  information),  the  INS  will  provide  open 
loop,  carrier  type  estimates  (derived  primarily  from  the  esti 
mates  of  the  aircraft's  velocity)  and  aircraft  dynamics  to 
the  code  loop.  The  INS  open  loop  information  results  in  con¬ 
siderable  improvement  in  the  second-order  code  loop's  track¬ 
ing  capability.  The  tracking  bandwidths  (of  both  code  and 
carrier  loops),  for  an  integrated  GPS/INS,  allow  substan¬ 
tially  higher  jamming  tolerance  and  improved  capability  to 
track  severe  aircraft  dynamics  as  compared  to  the  stand-alone 
GPS  UE  set  (5s 17-18) . 

UE  Satellite  Selection.  The  ultimate  objective  of  any 
satellite  selection  scheme  for  the  GPS  UE  is  to  pick  from  the 


available  satellites  the  set  which  minimizes  the  user’s  nav¬ 
igation  errors.  In  a  normal  environment  (without  jamming) 
the  optimum  set  of  four  satellites  to  track  is  determined  by 
the  geometry  between  the  user  and  the  in-view  satellites. 

For  example,  without  considering  an  unknown  time  bias,  the 
best  three  satellite  set  for  determining  the  user’s  position 
is  the  set  with  mutually  orthogonal  line-of -sight  (LOS)  unit 
vectors  with  the  origin  at  the  user's  position  extending  in 
the  direction  of  the  satellite.  The  orthogonal  geometry  pro¬ 
vides  maximum  information  along  all  three  coordinates  of  the 
user's  position.  If  the  satellites,  instead  of  being  aligned 
in  an  orthogonal  set,  are  grouped  or  clustered  closely  to¬ 
gether,  the  information  along  one  or  more  of  the  directions 
of  the  user's  position  is  degraded.  In  two  dimensions  for 
example,  the  user's  position  (the  "X”  or  aircraft)  in  the 
following  diagram  is  determined  by  the  intersection  of  the 
range  vector  to  each  satellite.  As  a  result  of  the  inaccu¬ 
racy  in  the  range  measurement  to  each  satellite,  the  inter¬ 
section  of  the  two  range  vectors  creates  an  area  of  uncer¬ 
tainty  within  which  the  user's  position  is  located  instead 
of  one  unique  solution  to  the  user's  position. 

The  shape  of  the  region  of  uncertainty  in  the  user's 
position  is  a  result  of  the  user -satellite  geometry.  For  a 
favorable  geometry  (when  satellite  LOS  vectors  are  at  or 
near  90°  separated  from  each  other)  the  region  of  user  posi¬ 
tion  uncertainty  is  approximately  a  square  with  dimensions 
defined  by  the  inaccuracies  of  the  measurements.  An 


unfavorable  geometry  creates  a  region  of  user  position  uncer¬ 
tainty  that  is  elongated,  indicating  increased  uncertainty 
in  user  position  along  the  direction  of  elongation.  Expand¬ 
ing  from  two  dimensions  to  three  dimensions  the  region  of 
user  position  uncertainty  changes  from  an  area  to  a  volume. 
Again,  a  favorable  satellite  LOS  geometry  will  place  the 
user  on  equal  distance  from  the  surface  of  the  uncertainty 
(or  error)  volume.  An  unfavorable  geometry  will  have  an 
error  volume  that  is  elongated  or  depressed  in  shape,  re¬ 
sulting  in  the  user's  position  having  greater  uncertainty 
in  the  direction  of  the  elongation  (the  direction  of  the 
grouped  or  clustered  satellites).  As  discussed  in  the  Space 
Segment  Section,  the  satellite  constellation  was  chosen  to 
minimize  the  clustering  or  grouping  of  satellites. 

The  time  bias  between  the  satellite  and  UE  sets'  clock 
adds  an  additional  range  error  to  each  satellite  measurement. 
If  the  time  bias  is  not  determined,  the  size  of  the  three 
dimensional  error  volume  would  grow  substantially.  To  deter¬ 
mine  and  remove  the  range  error  due  to  the  time  bias,  a 
fourth  satellite  range  measurement  is  needed.  Solving  for 
the  time  bias,  however,  adds  another  unfavorable  geometry 
condition.  The  unfavorable  time  bias  geometry  occurs  when 
all  four  satellites  are  close  to  or  on  the  surface  of  a  cone 
(7i4).  When  this  unfavorable  geometry  occurs  the  informa¬ 
tion  along  one  or  more  of  the  coordinates  of  the  user’s  posi¬ 
tion  becomes  redundant  with  the  time  bias,  thus  effectively 
reducing  the  number  of  range  measurements  to  three  and 


preventing  calculation  of  the  time  bias.  The  loss  of  observ¬ 
ability  of  the  time  bias  rapidly  increases  the  user's  posi¬ 
tion  uncertainty.  The  unfavorable  cone  geometry  is  a  func¬ 
tion  not  only  of  the  satellite  constellation,  but  the  specific 
user’s  location. 

The  UE  set  uses  a  figure  of  merit  known  as  geometrical 
dilution  of  precision  (GDOP)  to  select  from  the  available 
satellites  the  four  satellite  set  with  the  best  geometry  for 
tracking. 

GDOP  =  SQRT  ( TRACE ( ( H^H ) ~ 1 ) )  (7) 

The  visibility  matrix  (or  H  matrix)  is  formed  by  calculating 
the  direction  cosines  of  the  line-of -sight  unit  vector  from 
the  user  to  each  in-view  satellite  (expressed  in  user  local 
level  coordinates)  (6:144-145).  Each  row  of  the  H  matrix  is 
composed  of  the  three  direction  cosines  of  the  LOS  unit  vec¬ 
tor  and  the  fourth  element  for  all  satellites  is  a  one  to 
account  for  the  equal  uncertainty  in  range  due  to  satellite- 
user  clock  time  bias  for  each  satellite 

VD  V1’  AU(1)  1 

Ag(2)  A n(2)  Au(2)  1 

■  •  •  • 

H  =  (8) 

•  •  ♦  • 

•  •  *  • 

Vj)  Au(j)  1 


where 


A^*).  A^.),  A^(*)  are  the  respective  direction  co¬ 
sines  of  the  line-of-sight  unit  vector  to  each  satel¬ 
lite  expressed  in  user-local  level  (EAST-NORTH-UP)  co¬ 
ordinates 

j  =  the  number  of  in-view  and  available  satellites 

Due  to  variations  in  geographical  terrain,  all  satellites 
below  5°  elevation  from  the  horizon  are  deleted  from  possible 
satellite  selection  (the  5°  angle  is  known  as  the  mask  angle) 
(4:E9.3.2).  The  GDOP  figure  of  merit  (Equation  (7))  calcu¬ 
lates  a  scalar  value  for  each  possible  satellite  set.  The 
satellite  set  with  the  lowest  GDOP  is  the  best  user-satellite 
geometry  to  provide  the  smallest  navigation  errors  (7:4). 
Other  useful  figures  of  merit  are  derived  from  GDOP.  The 
position  dilution  of  precision  (PDOP)  is  calculated  using 
Equation  (7)  but  with  the  H  matrix  composed  of  the  direction 
cosines  to  the  satellites  only  (time/geometry  conditions  de¬ 
termined  by  the  fourth  column  of  l's  are  deleted).  PDOP  is 
used  extensively  in  satellite  constellation  design  and  analy¬ 
sis  without  having  to  specify  the  user's  position  and  trajec¬ 
tory.  Horizontal  dilution  of  precision  is  a  horizontal  plane 
accuracy  descriptor  like  CEP.  Vertical  dilution  of  preci¬ 
sion  is  valuable  in  terrain  avoidance.  The  next  chapter  will 
examine  GDOP  and  other  satellite  selection  criterion  in  more 
detail . 

UE  Adaptive  Antenna  Theory .  The  objective  behind  using 
an  adaptive  antenna  over  a  regular  antenna  is  to  limit  the 


noise  passed  on  to  the  satellite  tracking  loops,  and  thereby 
to  increase  the  performance  of  the  set.  This  capability  is 
especially  critical  in  a  jamming  environment.  The  adaptive 
antenna  designs  used  with  the  high  performance  GPS  UE  are 
either  a  beam  steering  or  a  null  steering  antenna.  Both  de¬ 
signs  are  phased  array  antenna  systems.  Characteristics  of 
a  phased  array  antenna  system  are  the  ability  to  provide 
rapid  inertialess  scanning  for  high  speed  angular  coverage, 
a  large  power  handling  capability,  and  their  suitability  to 
environments  where  actual  movement  of  the  antenna  is  diffi¬ 
cult  or  impossible  (9s 1).  Since  the  satellites  may  be 
tracked  anywhere  from  the  horizon  to  the  zenith,  both  antenna 
designs  must  provide  total  coverage  of  the  upper  hemisphere. 

The  problem  for  the  antenna  is  to  track  the  satellite  tra¬ 
jectories  while  limiting  jamming  power  (by  either  ground  or 
airborne  jammers)  being  passed  on  to  the  satellite  tracking 
loops . 

Beam  Steering  Antenna .  The  beam  steering  antenna 
creates  and  points  high  gain  beams  toward  each  satellite 
being  tracked,  providing  maximum  gain  in  the  direction  of 
the  satellite  and  very  low  gain  in  all  other  directions  (10:6). 
The  pointing  information  for  the  antenna  beams  is  derived 
initially  from  the  UE  estimate,  which  may  be  aided  with  INS 
information  (for  the  integrated  GPS/INS).  Once  tracking  of 
the  satellites  has  commenced,  pointing  directions  are  pro¬ 
vided  by  the  tracking  loops  and  navigation  information  (to 
compensate  for  vehicle  dynamics).  The  antenna  beam  width  and 


gain  determine  the  required  accuracy  of  the  pointing  informa¬ 
tion  and  the  jamming  tolerance  of  the  beam  steering  antenna. 
The  beam  steering  antenna  has  the  advantage  of  being  able  to 
overcome  moderate  jamming  levels  by  pointing  narrow,  high 
gain  beams  at  the  satellites  in  order  to  isolate  the  GPS 
signal  from  jamming  noise.  The  disadvantages  are  the  large, 
bulky  array  apertures  required  to  form  narrow  beams  at  the 
L  band  frequencies,  the  high  computational  requirements  for 
beam  steering  in  a  high  performance  aircraft,  and  the  unde¬ 
sirable  amplification  of  jamming  noise  in  the  direction  of 
the  satellite  (10s6). 

Null  Steering  Antenna .  The  null  steering  antenna  uses 
the  fact  that  the  satellite  signal  is  30  decibels  below  ther¬ 
mal  noise  (10:6).  Therefore,  any  signal  that  rises  above 
thermal  noise  must  be  interference.  The  antenna  algorithm 
acts  to  distribute  its  nulls  to  minimize  the  total  power  out 
of  the  antenna  (6:139).  This  process  is  accomplished  by 
combining  the  outputs  of  several  individual  antenna  elements 
through  phase  and  amplitude  controls  to  minimize  the  power 
out  of  the  antenna  (10:6).  The  advantages  of  the  null  steer¬ 
ing  antenna  are  its  smaller  array,  better  jammer  suppres¬ 
sion,  and  simplicity  of  operation  (compared  to  the  beam  steer 
ing  antenna).  The  disadv?  ges  of  the  null  steering  antenna 
are  that  the  number  of  nulls  it  can  create  is  limited  to  the 
number  of  antenna  array  elements  minus  one  (i.e.,  a  seven 
element  array  has  a  maximum  of  six  nulls)  and  that  in  an  ef¬ 
fort  to  null  out  strong  jammers,  other  "spurious”  nulls  may 


be  formed  that  null  a  satellite  signal  (10s6,  6:139-141). 

The  last  disadvantage  can  be  limited  by  constraining  the 
power  inversion  process  from  forming  a  null  in  the  direction 
of  a  satellite  (10:6,  6:140).  The  null  steering  antenna  will 
be  discussed  in  greater  detail  in  Chapter  IV,  and  in 
Appendix  B. 


Summary 

This  chapter  has  reviewed  the  three  segments  of  the 
GPS  Program  with  particular  emphasis  on  the  Satellite  and 
User  Segments.  The  Satellite  Segment  discussed  the  satellite 
constellation  and  the  GPS  satellite  transmitted  signal.  The 
User  Segment  presented  the  theory  and  operation  of  the  User 
Equipment.  The  User  Equipment  operation  included  a  discus¬ 
sion  of  the  GDOP  method  for  satellite  selection  and  adaptive 
beam  and  null  steering  antenna  designs  used  with  high  per¬ 
formance  GPS  User  Equipment.  The  next  chapter  will  further 
analyze  the  GPS  UE  satellite  selection  process. 


Introduction 


III  GPS  UE  Satellite  Selection 

The  standard  procedure  for  UE  satellite  selection  is 
based  solely  on  geometrical  dilution  of  precision  (GDOP)  (as 
discussed  in  the  User  Segment  Section  of  Chapter  II).  Nav¬ 
igation  measurements  are  dependent  on  three  factors  (exter¬ 
nal  to  the  UE  signal  processor):  the  geometry,  the  incoming 
signal,  and  the  noise  (11*1) .  GDOP  provides  a  figure  of 
merit  to  determine  the  best  geometry,  but  does  not  address 
the  incoming  t^yiial  and  noise.  An  optimum  satellite  selec¬ 
tion  process  must  use  every  available  piece  of  information 
available  to  determine  the  satellite  set  which  will  provide 
the  smallest  navigation  errors  over  a  particular  time  period. 
The  following  sections  of  this  chapter  will  discuss  this 
problem  and  possible  solutions  to  the  problem  in  much  greater 
detail.  The  next  section  will  review  previous  work  in  the 
area  of  GPS  UE  satellite  selection.  The  section  following 
the  Literature  Review  will  present  a  tutorial  on  GDOP  and 
improvements  and  extensions  to  GDOP  (from  a  more  rigorous 
stochastic  analysis) .  The  material  in  this  section  is  an 
excerpt  from  a  report  authored  by  Dr.  William  L.  Brogan, 
University  of  Nebraska-Lincoln,  under  contract  to  the  Avion¬ 
ics  Laboratory  of  the  Air  Force  Wright  Aeronautical  Labora¬ 
tories  (AFWAL)  at  Wr ight-Patterson  AFB,  Ohio  (11).  The 
work  in  this  section  coupled  with  the  development  in  the 
following  section  will  provide  the  reader  with  a  thorough 


understanding  of  the  concepts  and  theory  behind  the  UE  sat¬ 
ellite  selection  problem  and  possible  solutions  to  that  prob¬ 
lem. 


Literature  Review 

ADGINT  Study.  In  February  1981,  a  summary  of  the  work 
completed  by  the  Charles  Stark  Draper  Laboratory  on  the  "Ad¬ 
vanced  GPS /Inertial  Integration  Technology  Program  (ADGINT)", 
was  published  by  the  AFWAL  Avionics  Laboratory  (6).  In  that 
work  a  weighted  GDOP  satellite  selection  criteria  was  de¬ 
veloped  and  tested.  The  weighted  GDOP  criteria  was  a  sim¬ 
plification  of  the  a-posteriori  error  covariance  P^  equations 

Px  =  (Pq-1  +  HTR-1H)~1  (9) 


where  s 

H  is  the  user-satellite  geometry  matrix 
Pq  is  the  a  priori  error  covariance 
R  is  the  covariance  of  the  measurement  noise 

The  new  satellite  selection  criterion  will  select  the  best 
satellite  set  for  maximum  reduction  in  the  navigation  error 
covariance.  The  initial  error  covariance  P^  was  dropped  from 
equation  (9),  because  the  new  satellite  set  would  not  affect 
Pq,  thus  equation  (9)  reduces  to 

P.  =  (HTR_1H)“1  (10) 


T  —1  —1 

Therefore,  by  minimizing  TR (H  R  H)  the  a  posteriori  error 
covariance  is  minimized  (6s 145).  The  further  justifica¬ 
tion  of  dropping  from  equation  (9)  was  that  would  be 
very  large  (therefore  the  Pq-^  would  be  negligible)  and  that 
equation  (10)  would  result  in  a  smaller  number  of  computa¬ 
tions  than  equation  (9)  (6s 145). 

The  assumption  of  uncorrelated  pseudo-range  measure¬ 
ments  from  different  satellites  yields  a  diagonal  R  matrixs 


The  elements  of  R  (P^  P2,  P3,  and  P4)  are  functions  of  an 
estimate  of  the  received  signal -to-noise  ratio.  In  the 
ADGINT  study,  a  signal-to-noise  estimation  algorithm  was 
derived  from  outputs  of  a  null  steering  antenna  in  order  to 
aid  an  extended  range  receiver  tracking  design.  Therefore, 
the  signal-to-noise  estimate  was  readily  available  for  use 
in  the  satellite  selection  criteria  (6s29-33).  The  diagonal 
R  matrix  results  in  easy  matrix  inversion  using  a  very  small 
number  of  calculations. 


range  from  zero  to  one  depending  on 


The  weights,  W^,  in  R 
the  signal-to-noise  ratio.  The  upper  limit  of  W^  is  the 
signal-to-noise  level  at  which  further  performance  improve¬ 
ments  are  no  longer  advantageous  (compared  to  the  added  com¬ 
plexity  of  the  algorithm) .  The  lower  limit  of  W.  is  the 
signal-to-noise  level  at  which  acquisition  and  tracking  is 
no  longer  possible  (6s 146).  The  resultant  weighted-GDOP 
(WGDOP)  figure  of  merit  is, 

WGDOP  =  SQRT  (TRACE (HTR"1H)"1)  (13) 


The  study's  conclusions  of  the  performance  of  WGDOP  versus 
GDOP  in  a  GPS/INS  computer  simulation  determined  that  (6s 152): 

(1)  The  improved  signal-to-noise  power  of  the  satellite 
set  selection  by  the  WGDOP  criterion  more  than  compensates 
for  any  degradation  in  geometry. 

(2)  The  WGDOP  criterion  is  insensitive  to  small  changes 
in  signal-to-noise  power,  therefore,  only  a  crude  estimate 
of  signal-to-noise  power  is  needed  for  the  selected  set  to 
be  effective. 

(3)  The  WGDOP  criterion  should  be  used  more  frequently 
than  the  rate  used  by  GDOP  and  the  rate  would  be  dependent 
on  the  expected  jamming  threat. 

(4)  Most  of  the  benefits  of  WGDOP  are  derived  from 
eliminating  the  satellites  with  obviously  poor  signal-to-noise 
power.  A  selection  scheme  that  eliminated  satellites  with 
poor  signal-to-noise  power  and  then  used  the  standard  GDOP 
would  also  result  in  performance  improvements. 
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The  ADGINT  study  based  the  WGDOP  versus  GDOP  conclusions 
on  a  computer  simulation  with  one  particular  mission  in  which 
the  jamming  threat  was  encountered  at  the  end  of  the  mission. 
Encountering  the  jamming  threat  after  a  long  period  of  normal 
operation  would  result  in  a  small  a  priori  error  covariance 
PQ,  which  is  in  conflict  with  the  assumption  made  to  elimi¬ 
nate  P  from  the  satellite  selection  criteria.  Therefore, 

— o 

minimizing  WGDOP  may  not  in  fact  result  in  the  smallest  a- 
posteriori  error  covariance  P^ . 

In  March  1981,  a  short  study  on  GPS  UE  satellite  selec¬ 
tion  was  published  by  the  Reference  Systems  Branch  of  the 
AFWAL  Avionics  Laboratory  (7).  In  that  study  a  combined 
jamming-to-signal  (J/S)  estimate  and  GDOP  criteria  was  used 
for  satellite  selection.  The  study  sought  to  prove  the  value 
of  incorporating  a  generalized  form  of  J/S  estimate  into  the 
satellite  selection  criteria,  and  for  this  reason  a  generic 
UE  receiver  and  adaptive  antenna  were  simulated  (the  generic 
UE  receiver  and  adaptive  antenna  were  used  to  make  resulting 
performance  improvements  applicable  to  all  GPS  UE  receiver 
and  adaptive  antenna  designs  regardless  of  the  methods  for 
processing  the  satellite  signals  inside  the  set) .  Naviga¬ 
tion  performance  was  judged  to  be  improved  if  the  best  GDOP 
set  included  one  or  more  satellites  above  a  J/S  threshold  of 
70  decibels  while  an  alternate  satellite  set  existed  that 
had  all  satellites  below  the  70  db  J/S  level  and  had  a  GDOP 
value  below  5.0.  This  procedure  for  satellite  selection  is 
very  similar  to  the  procedure  proposed  in  conclusion  (4)  of 
the  preceding  section.  The  study  tested  the  combined  J/S 


and  GDOP  versus  GDOP-only  satellite  selection  criterion 
against  a  stylized  area  of  jammers  (ten  rows  of  ground  jam¬ 
mers,  paralleling  the  forward  edge  of  the  battle  area  and 
extending  80  kilometers  deep  into  enemy  territory  and  stretched 
beyond  the  aircraft  line-of-sight  in  both  directions  ( 7 1  7 ) ) 
broadcasting  one  kilowatt  each.  Five  different  aircraft  mis¬ 
sions  into  enemy  territory  were  simulated  (Close  Air  Support, 
Escort/Intercept,  Attack  Helicopter,  Preplanned  Deep  Strike, 
and  Intratheater  STOL  Transport) .  The  results  of  the  simula¬ 
tion  for  a  24  satellite  constillation  showed  an  average  per¬ 
formance  improvement  across  all  5  aircraft  missions  of  42% 
from  a  10°  mask  angle  and  15%  for  a  20°  satellite  mask  angle 
(7s 15).  The  results  for  an  18  satellite  constellation  showed 
an  average  performance  improvement  across  all  5  aircraft  mis- 

O  o 

sions  of  15%  for  a  10  mask  angle  and  3%  for  a  20  mask 
angle  (7s 16).  Take  note  here  that  the  18  satellite  constel¬ 
lation  referred  to  is  a  predecessor  of  the  present  18  satel¬ 
lite  constellation.  The  difference  between  the  two  constel¬ 
lations  is  that  originally  six  satellites  would  be  unevenly 
spaced  within  one  orbital  plane,  with  a  total  of  three 
orbital  planes,  and  zero  phasing  of  satellites  from  one 
plane  to  another  (phasing  is  discussed  in  Chapter  II)  (4:E9.3.1 
E9.3.3).  A  major  concern  identified  by  this  study  is  the  im¬ 
portance  of  satellite  availability  at  the  satellite  selection 
time.  The  24  satellite  constellation  provided  an  average  of 
8  visible  satellites  with  70  alternative  satellite  sets  pos¬ 
sible.  The  18  satellite  constellation  provides  5-6  visible 


satellites;  6  satellites  give  15  alternative  satellite  sets 
while  5  satellites  result  in  only  5  alternative  satellite 
sets.  The  satellite  visibility  effect  has  some  direct  con¬ 
sequences  on  any  alternate  satellite  selection  scheme: 

(1)  The  added  complexity  and  computations  of  an  alter¬ 
nate  satellite  selection  scheme  may  not  be  justified  if  there 
are  few  or  no  alternative  satellite  sets. 

(2)  The  number  of  computations  for  a  complex  satellite 
selection  scheme  would  also  be  greatly  reduced  with  the  num¬ 
ber  of  visible  satellites.  The  reduction  in  the  number  of 
computations  would  allow  a  satellite  selection  scheme  with 
greater  complexity  and  larger  performance  improvement. 

Clearly  (1)  and  (2)  pose  a  degree  of  conflict  in  design  of 
an  alternate  satellite  selection  criteria  that  must  be  ad¬ 
dressed  as  well  as  any  performance  improvements. 

Satellite  Visibility.  The  present  18  satellite  con¬ 
stellation  replaced  the  old  satellite  constellation  because 
of  increased  performance  of  the  uniform  3  satellites  in  6 
different  orbits  constellation  over  the  non-uniform  6  sat¬ 
ellites  in  3  different  orbits  (See  Reference  4:E9.3.1-E9.3.8) 
For  further  analysis  of  the  satellite  visibility  of  the  pre¬ 
sent  18  satellite  constellation  the  software  used  in  the 
AFWAL  Avionics  Laboratory  study  was  obtained  and  modified  to 
provide  Figure  9  (Ref.  7). 

For  a  further  discussion  of  the  satellite  visibility 
software  used  to  create  the  information  for  Figure  9,  see 
Appendix  A.  The  results  of  the  satellite  visibility  study 


using  the  present  18  satellite  constellation  show  that  for  a 
mask  angle  of  5°  (which  is  typically  used  in  the  GPS  UE)  an 
average  of  6  satellites  would  be  available  for  tracking  over 
a  24  hour  period. 


TABLE  I 


Satellite  Availability 

for  a  24  Hour  Period 

Mask  Angle 

Average  Number  of 
Visible  Satellites 

0° 

7.0 

5° 

6.0 

10° 

5.6 

15° 

4.9 

20° 

4.2 

25° 

4.0 

These  results  are  in  agreement  with  the  "NAVSTAR  Global 
Positioning  System  Six-Plane  18-Satellite  Constellation," 
data  (4:E9.3.1-E9.3.8) . 

AFWAL  Brogan  Study.  In  September  1981,  the  AFWAL 
Avionics  Laboratory  published  the  study  by  Dr.  William  Brogan 
that  included  an  optimal  stochastic  approach  to  an  alternate 
satellite  selection  criterion.  Within  this  study,  a  complete 
stochastic  analysis  of  the  GDOP  concept  for  selecting  satel¬ 
lites  for  navigation  measurements  was  presented.  In  that  de¬ 
velopment  of  an  alternate  satellite  selection  criterion,  not 
only  was  the  satellite  visibility,  user-satellite  geometry. 


and  jamming  noise  at  the  receiver  considered,  but  also  the 
a  priori  error  covariance  P  that  was  deleted  from  the  ADGINT 


weighted-GDOP  criteria.  The  value  of  to  determine  direc- 
tion(s)  of  the  maximum  "need  for  information"  in  the  satel¬ 
lite  selection  process  adds  a  great  deal  o f  insight  into 
which  of  the  possible  satellite  range  measurements  can  result 
in  the  smallest  overall  navigation  errors.  Before  discussing 
the  alternate  satellite  selection  criterion  that  will  be 
implemented  and  tested  in  this  report,  a  more  detailed  re¬ 
view  of  Dr.  Brogan’s  study  is  essential.  Due  to  the  clarity 
and  insight  of  Dr.  Brogan's  work  the  following  section  will 
present  excerpts  from  Dr.  Brogan's  study  (11). 

The  Cost  Criterion  for  UE  Satellite  Selection  (11) 

An  Introduction  to  GDOP  Using  Least-Squares  (11s  3-7 ) . 
Consider  a  vector  of  measurement  components,  Z,  corrupted  by 
additive  noise  N  and  linearly  related  to  the  state  vector  X. 

Z  =  H  X  +  N  (14) 

.  A 

An  estimate  of  X,  called  X,  is  sought  from  the  given  measure¬ 
ment  vector  Z.  One  approach  to  this  problem  is  the  least- 

a  .  A  ... 

squares  estimator  X.  That  is,  select  X  so  as  to  minimize 

the  sum  of  the  squares  of  the  components  of  the  difference 
A  a 

Z  -  HX.  This  X  would  then  minimize 

Jx  =  (Z-HX)T(Z-HX)  (15) 

A 

Differentiating  with  respect  to  X  and  setting  the  result 
to  zero  gives 
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(16 


-H~Z  +  H~HX  =  0 


so  that  the  well  known  least-squares  estimate  is 


A  T  — 1  T 

X  =  (HH)  UZ 


The  quality  of  this  estimate  can  be  assessed  by  forming  the 
error 


AJ  A 
X  =  X  -  X 


Using  equation  (14)  and  (17)  gives 


X  =  (hth)~1hthx  +  (HTH)"1HTN  -  X 


=  (hth)_1htn 


If  the  noise  N  has  zero  mean,  taking  the  expected  value 
E[  ]  of  equation  (19)  shows  that  the  estimator  error  X 
also  has  a  zero  mean.  The  covariance  of  the  estimation 
error  can  be  related  to  the  covariance  of  the  noise  N  bys 


Ttt\  “IttT—  »,%tT  tt/t»Ttt\—  1 


E[ja  ]  =  (HH)  P  NN  H(HH) 


If  all  components  of  N  are  pairwise  uncorrelated  and  have 


unit  variance,  i.e.. 


eUn.N  .]  =  [6  .  .]  = 

3  3  0  if  i  /  j 

then  E[NNT]  =  I  so  that 

e[xx]t  =  (HTH)“1(HTH)(HTH)“1 

=  (HTH)~1  (21) 

T  —1 

The  matrix  (H  H)  is  called  the  GDOP  matrix.  Various  scalar 
figures  of  merit  can  be  derived  from  this  matrix,  and  are 
frequently  used.  For  example,  the  trace  of  (H  H)  ,  is  the 
sum  of  the  squares  of  all  error  components  when  E[N.N  .]  »  [6-  0 

l  ^  J  ^  J 

The  square  root,  [TRACE  (HTH)“13'S,  is  usually  called  GDOP. 

The  square  root  of  the  sum  of  all  diagonal  terms  which  cor¬ 
respond  to  position  components  is  called  PDOP.  If  the  hor- 
izontal  plane  diagonal  components  of  (H  H)  are  summed,  the 

square  root  is  called  HDOP.  In  this  study,  major  attention 

T  -1 

is  directed  to  TRACE (H  H)  .  0 mission  of  the  square  root  will 

not  alter  any  conclusions  regarding  the  minimization  of  GDOP. 

It  is  clear  that  all  GDOP-related  performance  measures 
indicate  the  error  in  an  estimated  navigation  quantity  "per 
unit  of  measurement  noise"  covariance.  Deviations  from  this 
assumption  will  be  considered  in  the  next  subsection. 

All  of  the  above  GDOP-related  measures  depend  solely 
on  the  geometry  matrix  H.  Smaller  GDOP  values  indicate 
stronger  or  more  robust  geometric  solutions  to  the  estimation 


problem.  For  these  reasons,  when  some  freedom  exists  in 
the  choice  of  measurements,  good  (i.e.,  small)  GDOP  is  often 
used  as  the  selection  criteria.  It  will  be  shown  later  that 
this  criterion  is  not  always  the  best  choice. 

Weighted  Least  Squares  Solution  (11:5-7).  The  "per 
unit  noise”  concept  inherent  in  GDOP  is  not  so  useful  (or 
more  bluntly,  invalid)  when  certain  measurement  components 
are  much  noisier  than  others.  When  a  choice  exists  between 
two  possible  sets  of  measurements,  the  set  with  poorer  GDOP 
may  be  preferable  if  they  are  of  sufficiently  lower  measure¬ 
ment  noise.  This  trade-off  was  pointed  out  in  conjunction 
with  an  Integrated  GPS/Inertial  Navigation  System  (12). 

There,  the  desire  to  achieve  good  GDOP  suggested  selection 
of  range  measurements  separated  by  angles  near  90°.  However, 
this  choice  caused  an  increase  in  atmospheric  diffraction 
errors  in  those  measurements  with  low  elevation  angles . 
Another  potentially  more  severe  problem  is  non-uniform  noise 
in  a  jamming  environment.  Some  of  the  highly  directional 
antennas  may  receive  large  amounts  of  noise  pollution.  This 
degradation  may  well  overcome  any  GDOP  advantages. 

If  the  measurement  noise  covariance  is  not  just  the 

unit  matrix,  but  rather  is  E[NNT]  -  I  then  the  GDOP  matrix 
.  2 

just  increases  by  the  scalar  a  .  If  all  potential  measure- 

.  2 

ments  have  this  same  variance  o  ,  the  choice  of  the  best 

measurement  set  will  still  be  the  best  GDOP  set.  The  multi- 

2  . 

plication  scalar  a  will  not  affect  the  relative  rankings. 

The  measurement  noise  is  called  non-uniform  when  dif¬ 
ferent  measurements  have  different  noise  levels,  as  indicated 


by  their  variances.  In  this  case  a  weighted  least-squares 
approach  to  estimation  is  often  used  (19).  The  quadratic 
form  of  equation  (15)  is  modified  by  inserting  a  weighting 
matrix  W. 

J2  =  (Z-HX)TW(Z-HX)  (22) 

It  is  customary  to  select  the  weighting  matrix  as  the  inverse 
of  the  noise  covariance  matrix  (i.e.,  W  =  R-1  where 
R  =  E[NNT]).  This  choice  weights  the  accurate  measure¬ 
ments  more,  the  noisy  ones  less.  Minimizing  this  modified 
cost  function  leads  to  a  modification  of  equation  (15). 

X  =  (hV1H)'1hV1Z  (23) 

and  a  corresponding  modification  of  equation  (20)  or  (21) 

E[XXT]  =  (HTR-1H)~1  (24) 

Equation  (24),  and  generalizations  of  it,  are  Known  as  noise- 
weighted  GDOP .  It  is  presented  here  to  contrast  against  the 
traditional  GDOP,  but  will  be  discussed  further  in  the  next 
subsection. 

Thus  far,  no  mention  has  been  made  of  the  knowledge 
about  the  navigation  states  which  may  exist  prior  to  making 
GPS  measurements.  In  almost  all  real  situations  prior  know¬ 
ledge  of  the  navigation  states  will  exist  because  of  previous 


measurements  or  for  other  reasons.  The  uncertainty  volume 
of  X  need  not  be  spherical  (an  uncertainty  volume  of  X  is  a 
set  of  points  enclosed  within  an  equi-probability  surface) . 

An  example  of  an  error  volume  would  be  the  case  where  X,  Y, 
and  Z  components  of  position  uncertainty  are  equal.  In  two 
dimensions,  since  X  and  Y  are  equal  the  conic  is  a  circle 
with  a  radius  equal  to  X  and  Y.  Expanding  to  three  dimen¬ 
sions,  the  equi-probability  ellipsoid  is  a  sphere  with  a 
radius  equal  to  X,  Y,  and  Z.  Some  components  of  X  can  have 
much  more  uncertainty  associated  with  them,  leading  to  elon¬ 
gated  ellipsoidal  uncertainty  volumes.  Those  components 
with  the  largest  uncertainty  are  said  to  have  "maximum  need" 
for  improvement.  A  GPS  measurement  which  provides  maximum 
information  along  the  elongated  axis  will  likely  lower  nav¬ 
igation  uncertainty  the  most  regardless  of  classical  GDOP 
considerations . 

The  treatment  of  a  priori  information  in  estimation  can 
be  approached  in  several  ways.  One  way  is  to  modify  equation 
(15),  this  time  to 

J3  =  (z-recrR  ^Z-ffiC)  +  (Xq-X)  Pq  (Xq-X)  (25) 

where  5^  represents  an  a  priori  estimate  of  X  and  the  covari¬ 
ance  of  5^  is  (13).  Minimizing  amounts  to  a  trade-off 
between  measurement  residuals  and  deviations  from  X Q.  This 
approach  leads  to  a  modified  estimator 


(26a) 


x  =  ChTR  "  1h  +  p0"1]“1[hV1z  +  Pq’1^] 

=  Xq  +  PqHT[HP_HT  +  R]_1[Z  -  HXp~i  (26b) 

Form  (26a)  is  most  commonly  used  in  weighted  least-squares 
analysis,  while  form  (26b)  is  one  of  the  well  Known  Kalman 
Filter  equations  (Ref.  14).  The  corresponding  forms  of  the 
a  priori  estimation  error  covariance  matrix  are 

/v  a/ 

E[XXT]  =  P  =  [HTR"1U  +  Pq-1]-1  (27a) 

=  P0-P0HT[HPoHT  +  R]-1^  (27b) 

Note  that  equation  (26a)  looks  most  like  the  GDOP 

matrix  in  that  it  reduces  to  equation  (24)  if  P^  approaches 

infinity  (very  large  navigation  position  uncertainties)  and 

it  reduces  to  equation  (21)  if,  in  addition  to  Pq  being  very 

2 

large,  R  =  I  or  R  =  la 

Generalizations  and  Extensions  of  GDOP  (Ilil5-16).  A 
major  part  of  the  Brogan  study  was  devoted  to  the  question 
of  how  best  to  select  measurements  at  a  given  point  in  time. 
This  static  measurement  selection  problem,  so  called  because 
only  one  point  in  time  is  considered,  is  a  useful  simplifica¬ 
tion  to  the  dynamic  measurement  selection  problem  for  theo¬ 
retical  development  of  extensions  to  GDOP. 

In  the  GPS  problem  we  are  concerned  with  three  posi¬ 
tion  components  and  a  clock  bias,  which  converts  to  an 
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equivalent  user  range  error.  So  for  the  most  part  a  four 
component  state  vector  will  be  of  ultimate  interest.  Making 
the  usual  assumptions  about  mean  values  being  zero,  atten¬ 
tion  turns  to  the  second  statistical  moments  of  the  estima¬ 
tion  errors,  i.e.,  the  covariance  terms.  Subject  to  the 
usual  assumptions ,  equation  (26a)  or  equation  (26b)  tells 
the  entire  story.  Given  an  a  priori  state  geometry  matrix 
H  and  a  noise  covariance  matrix  R,  the  a  priori  covariance 

"  A 

is  uniquely  computable. 

It  is  awkward  to  us%  a  4  x  4  matrix  to  describe  the  ac¬ 
curacy  or  effectiveness  of  a  given  set  of  measurements.  An 
attractive  alternative  is  to  use  a  single  scalar  figure  of 
merit  to  describe  the  result.  In  the  past,  various  terms 
such  as  circular  error  probable  (CEP)  or  spherical  error 
probable  (SEP)  have  been  used.  The  a  posteriori  covariance 
matrix  P^  can  be  viewed  geometrically  as  a  hyperellipsoid 
(similar  to  the  error  volume  discussed  earlier  but  expanded 
to  a  larger  dimension)  with  principal  axes  oriented  in  space 
according  to  the  directions  of  the  orthonormal  eigenvectors 
— i  of  Pi  and  t^ie  dimensions  of  the  semimajor  axes  as 

the  eigenvalues  of  P^.  The  interior  of  this  ellipsoid 
can  be  thought  of  as  a  one-sigma  error  volume,  and  in  this 
respect  it  is  related  to  CEP  and  SEP  although  much  more 
complicated.  Two  very  simple  scalar  measures  derived  from 
P^  can  be  used.  The  determinant  of  P^,  is  equal  to  the 
product  of  the  eigenvalues  (15) 


and  is  therefore  related  to  the  volume  of  the  elipsoid.  This 
measure  is  not  preferred  here  for  two  reasons.  First,  deter¬ 
minants  are  not  as  easy  to  evaluate  as  the  trace,  which  is 
used  in  the  next  performance  measure.  Second,  and  more 
important,  det(P^)  =  0  does  not  mean  the  error  is  zero. 

All  it  means  is  that  one  or  more  dimensions  of  the  hyper¬ 
ellipsoid  has  degenerated  to  zero,  but  very  large  errors  in 
other  coordinates  may  still  exist.  The  det(P1)  is  not  very 
descriptive  or  discriminating  in  this  regard. 

The  preferred  scalar  performance  measure  is  based  on 


the  trace  of  P^ 


N 


TR(P-)  =  E  P. 
1  i=l  1 


ll 


N 


=  E  A. . 


i=l  1 


(29) 


The  trace  is  simply  computed  as  the  Siam  of  the  diagonal  ele¬ 
ments,  but  it  is  known  that  it  also  equals  the  sum  of  the 
eigenvalues  (15).  Therefore  TR(P^)  can  not  equal  zero  unless 
all  error  variance  components  are  zero.  Also  the  trace  of 
P^  is  analogous  to  the  radius  squared  value  of  a  represen¬ 
tative  one  sigma  hyper-sphere  whereas  det(P1)  is  the  volume 
of  the  hyper-sphere.  The  trace  functions  just  described  are 
natural  extensions  of  GDOP.  In  the  following  section  two 
extensions  of  GDOP  will  be  considered:  noise-weighted  GDOP 
and  the  Cost  Criterion. 


The  Cost  Criterion  for  Alternate  Satellite  Selection. 
In  the  second  section  of  this  chapter,  equation  (24)  was 
called  the  noise-weighted  GDOP.  The  trace  of  equation  (24) 
may  also  be  called  noise-weighted  GDOP. 


rT„-ltI-i-l 


Noise-Weighted  GDOP  =  TRACE  [h R  h] 


(30) 


Noise-Weighted  GDOP  penalizes  noisy  measurements  in  the  satel 
lite  selection  process.  As  stated  in  the  first  section  of 
this  chapter,  in  the  discussion  of  WGDOP,  Noise-weighted 
GDOP  assumes  that  can  be  neglected  from  equation  (30) 

without  loss  of  performance.  For  Pq-1  to  be  neglected  would 
mean  that  P^  is  very  large,  which  could  be  expected  very 
early  in  a  mission  (before  the  accuracy  of  GPS  has  been 
achieved)  or  after  a  long  period  of  GPS  unavailability 
(which  could  occur  from  extended  periods  of  severe  jamming) . 
For  most  examples  of  interest  this  simply  will  not  be  the 
case  (11:18).  Therefore,  Pq-*  should  be  included  in  an  opti¬ 
mum  satellite  selection  criteria,  as  in  equation  (27a). 
Equation  (27a)  was  one  of  two  forms  of  the  a-posteriori 
estimation  error  covariance  matrix.  The  trace  of  equation 
(27a)  will  henceforth  be  known  as  the  cost  criterion. 


Cost  Criterion  =  TRACE (P^ ) 


=  TRACE  [HTR-1H  +  P-1]"1  (31) 


Tv  \  v  v  \  V 


The  term  “cost  criterion"  is  chosen  to  reflect  the  fact  that 


this  equation  is  a  result  of  the  cost  function  J^,  equation 
(25).  The  cost  criterion  is  a  tradeoff  between  noise¬ 
weighting  of  the  user-satellite  geometry  with  the  directions 
of  the  "need  for  information"  contained  in  P^.  The  added 
computations  from  the  noise-weighted  GDOP  criterion  are  min¬ 
imal  because  should  already  be  known  by  the  UE  set. 

As  in  the  GDOP  error  volume  discussion  (in  the  UE  Sat¬ 
ellite  Selection  Section  of  Chapter  II)  as  the  GDOP  value 
increases  the  error  volume  becomes  elongated  and  distorts 
the  UE  sets'  navigation  accuracy.  A  similar  result  occurs 
by  including  as  part  of  equation  (31).  As  the  estimation 
error  covariance  hypervolume  elongates  P^  will  become  more 
dominant  in  determining  the  best  satellite  measurements. 

Then  as  the  estimation  error  covariance  hypervolume  equalizes 
the  (H  R  H)  term  would  become  more  dominant  in  satellite 
measurement  selection  (an  equalized  hyper-volume  would  be 
similar  to  a  sphere  in  three  dimensions). 

The  cost  criterion  is  a  sound,  thorough,  and  promising 
approach  for  an  optimal  solution  to  the  static  GPS  UE  satel¬ 
lite  selection  problem.  It  is,  however,  unclear  whether  the 
cost  criterion  will  result  in  the  smallest  mean  square  user 
position  error  for  the  dynamic  measurement  selection  prob¬ 
lem.  For  this  reason  the  cost  criterion  will  be  tested 
against  GDOP-only  criteria  by  means  of  an  integrated  GPS/INS 
computer  simulation  for  a  realistic  high  dynamic  (both  air¬ 
craft  and  jamming  dynamics)  flight  trajectory. 
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One  caveat  on  the  cost  criterion  is  that  it  is  one  ap¬ 
proach  for  an  optimal  satellite  selection  process,  not  the 
optimal  satellite  selection  process.  The  cost  criterion  was 


based  on  the  cost  function  (equation  (25)).  Additional 
numbers,  locations,  power  levels,  and  types  of  jamming  sig¬ 
nal  and  a  penalty  function  for  acquiring  satellites  with 
moderate  to  high  signal-to-noise  levels  could  be  added  to 
J^.  A  cost  function  with  this  type  of  additional  informa¬ 
tion  would  mean  a  substantial  increase  in  complexity  and 
computational  loading  that  may  be  unwarranted  for  the  limited 
satellite  visibility  of  the  18  satellite  constellation. 
Analysis  of  the  cost  criterion  versus  GDOP  performance 
should  indicate  whether  added  complexity  and  computational 
loading  of  an  alternate  satellite  selection  criterion  is 
warranted . 

Summary 

This  chapter  has  presented  a  detailed  analysis  of  the 
static  measurement  selection  problem  with  particular  emphasis 
on  GDOP  and  extensions  to  GDOP  for  satellite  selection.  One 
extension  to  GDOP  uses  estimates  of  the  measurement  noise 
to  each  satellite  and  the  a  priori  user  position  error 
covariance  matrix  as  well  as  satellite  geometry  information 
to  select  the  best  satellite  set.  This  extension  to  GDOP 
is  called  the  cost  criterion.  Analysis  of  the  cost  criterion 
in  a  static  measurement  selection  problem  has  indicated 
significant  improvement  in  navigation  accuracy  (11:36-37). 


It  is  unclear,  however,  whether  the  cost  criterion  will  re¬ 
sult  in  the  smallest  mean  square  user  position  error  for  the 
dynamic  measurement  selection  problem.  This  report  will 
implement  the  cost  criterion  in  an  integrated  GPS/INS  com¬ 
puter  simulation  to  analyze  the  cost  criterion  versus  GDOP 
navigation  performance  and  verify  the  results  of  the  static 
measurement  case. 


Overview 


The  Air  Force  Avionics  Laboratory  (AFAL)  developed  a 
digital  simulation  computer  program  for  an  integrated  GPS/ 

INS  aircraft  navigation  system  (IGI  simulator)  in  1976  (16s 
841).  The  IGI  simulator  was  a  result  of  a  number  of  in-house 
(AFWAL)  computer  simulations  in  support  of  the  Generalized 
Development  Model  (GDM)  of  GPS  User  Equipment.  The  GDM 
was  designed  to  be  used  as  a  flight  test  bed  to  evaluate 
high  anti- jam  system  techniques  for  military  applications 
and  to  expand  the  technology  base  for  GPS  User  Equipment. 

The  IGI  simulator  is  based  on  the  GDM  and  is  designed  to  sup¬ 
port  a  variety  of  GPS  User  Equipment  design  trade-off  studies 
and  mission  analyses  (16:841).  Thus  the  IGI  simulator  is  a 
logical  choice  for  evaluation  of  a  cost  criterion  versus 
GDOP  satellite  selection  criterion.  The  following  section 
of  this  chapter  will  analyze  the  functions,  math  models,  and 
operation  of  the  IGI  simulator.  Also,  the  modifications  to 
the  IGI  simulator  necessary  to  evaluate  the  cost  criterion 
versus  GDOP  satellite  selection  criterion  in  a  jamming  envi¬ 
ronment  will  be  presented  in  the  final  section  of  this  chap¬ 
ter.  The  level  of  detail  presented  in  this  chapter  will  be 
limited  to  the  necessary  material  for  an  understanding  of 
the  IGI  simulator’s  operation  and  performance  outputs  that 
will  be  analyzed  in  the  next  chapter. 


IGI  Description  (17s 3-5).  The  IGI  simulator  uses  a 
direct  linear  simulation  approach  in  which  the  differential 
equations  of  the  INS  algorithm  are  linearized  about  a  pre¬ 
computed  flight  trajectory  to  describe  the  propagation  of 
INS  errors  in  position,  velocity,  and  attitude.  The  INS 
indicated  data  is  then  processed  together  with  GPS  pseudo¬ 
range  and  delta  pseudo-range  measurements  in  a  linearized 
Kalman  filter  (21:23-67)  to  obtain  estimates  of  the  errors 
in  the  INS-indicated  position  and  velocity,  and  of  other 
variables  describing  error  sources  that  might  affect  naviga¬ 
tion  performance.  The  user’s  new  navigation  parameters  (the 
INS-indicated  data)  are  updated  by  combining  the  Kalman 
filter  error  estimates  with  the  old  INS-indicated  data. 

This  direct  linear  simulation  approach  was  chosen  as  a  com¬ 
promise  between  the  less  realistic  covariance  analysis  ap¬ 
proach  (where  the  equations  describing  all  of  the  state 
variables  are  a  linear  -combination  of  the  state  variables) 
and  the  more  computationally  demanding  non-linear  simulation 
approach  ( 16 : 842 ) . 

A  Functional  Look  at  IGI .  The  IGI  simulator  consists 
of  a  three-program  package;  PROFGEN,  IGI,  and  PLOTTER. 
PROFGEN  (Ref.  24)  is  an  aircraft  profile  generator  which 
generates  a  realistic  flight  trajectory  data  file.  This 
flight  trajectory  data  file  along  with  a  set  of  user-deter¬ 
mined  configuration  and  initialization  parameters  are  input 
to  IGI.  The  PLOTTER  program  processes  the  IGI  plot  output 


file  to  create  time  history  plots  of  navigation  errors.  The 
IGI  program  consists  of  the  main  code  (or  executive)  and 
three  levels  of  subroutines.  Level  one  subroutines  perform 
the  major  functions  in  IGI.  Level  two  and  three  subroutines 
provide  specialized  processing  of  data  (such  as  storing  a 
block  of  data  for  printed  or  plotted  output)  and/or  utility 
routines  (such  as  matrix  multiplication,  addition,  and  in¬ 
version)  .  The  executive  determines  when  and  which  of  the 
major  functions  must  be  performed.  The  major  functions  of 
IGI  are, 

(1)  Initialization  and  Configuration 

(2)  Satellite  Propagation  and  Selection 

(3)  State  (Truth  and  Filter)  Propagation 

(4)  Measurement  Update 

(5)  Store  Output  Data 

A  graphical  view  of  the  IGI  simulator  (input  data,  flight 
trajectory  file,  the  IGI  executive,  and  the  five  major  func¬ 
tions  within  the  IGI  program)  is  presented  in  Figure  10. 

The  initialization  function  reads  data  from  the  input 
data  file  to  determine  the  simulation  user’s  desired  config¬ 
uration  and  initialization  parameters.  For  instance,  the  in¬ 
put  data  file  specifies  the  number  of  truth  and  filter  states 
the  time  interval  for  storing  printed  and  plotted  output, 
the  time  interval  for  satellite  selection,  the  number  and 
type  of  external  measurements  (i.e.,  pseudo-range  and/or 
delta  pseudo-range  GPS  measurtments ) ,  and  the  statistics  of 
the  simulation’s  error  models. 
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Figure  10.  Structure  of  the  IGI  Simulation 


The  second  major  IGI  function  simulates  the  GPS  sat¬ 
ellite's  positions  and  velocities  at  the  simulation  time. 

The  satellite  selection  process  uses  the  GDOP  criterion  to 
determine  the  best  four  satellite  set  for  tracking  (See  UE 
Satellite  Selection  Section  of  Chapter  II).  Satellite  selec¬ 
tion  is  carried  out  once  every  five  minutes  to  insure  the 
best  satellite  set  is  being  tracked  (modeling  of  the  satel¬ 
lites'  positions  and  velocities  will  be  discussed  further 
in  Chapter  IV ) . 

The  state  propagation  function  for  both  the  truth  and 
filter  states  is  accomplished  by  continuously  updating  the 
truth  and  filter  states  with  INS-derived  data  alone.  State 
propagation  occurs  between  GPS  measurements  or  during  periods 
of  GPS  nonavailability  (as  would  be  expected  during  severe 
jamming  conditions) .  (A  more  detailed  discussion  will  be 
presented  in  the  IGI  Operation  Section  of  Chapter  IV). 

The  measurement  update  function  creates  the  true  measure 
ments  (satellite-user  ranges  and  range-rates)  then  corrupts 
the  true  measurements  with  various  error  sources  (such  as 
satellite-user  clock  bias,  atmospheric  delay,  ionospheric 
errors,  and  environmental  noise  errors)  to  create  the  four 
GPS  pseudo-range  and  delta  psdueo-range  measurements.  The 
measurement  information  is  used  to  correct  the  filter  state 
estimates  by  sequentially  using  Carlson's  upper  triangular 
square  root  algorithm  for  scalar  measurements  (19)  (discus¬ 
sed  in  more  detail  in  the  Satellite  Measurement  Generation 
and  Incorporation  Subsection  of  the  IGI  Operation  Section 
of  Chapter  IV ) . 


The  simulator  outputs  are  of  two  forms,  printed  and 
plot  data  output.  The  printed  output  echoes  the  data  from 
the  input  files  (both  the  flight  trajectory  and  configuration 
initialization  data  files)  and  the  program-generated  outputs 
before  and  after  a  measurement  update.  The  program-generated 
outputs  include  the  following  (17:68), 

(1)  True,  INS -indicated,  and  filter-estimated  whole 
value  user  position  and  velocity  variables. 

(2)  The  square  root  of  the  filter  covariance  diagonal 
elements . 

(3)  The  elements  of  the  filter  error  state  vector. 

(4)  The  elements  of  the  filter  estimation  error. 

(5)  The  elements  of  the  true  error  state  vector. 

(6)  The  satellites  in  view  at  the  satellite  selection 

time. 

(7)  The  set  of  four  satellites  with  the  minimum  GDOP. 

(8)  The  GDOP  value  of  the  optimum  satellite  set. 

The  simulation  plot  data  file  contains  a  series  of 

navigation  data  records.  A  navigation  data  record  contains 
the  time,  the  elements  of  the  filter  estimation  errors,  the 
elements  of  the  square  root  of  the  covariance  diagonal,  and 
the  INS-indicated  position  and  velocity  errors. 

IGI  Truth  and  Filter  States  (16:842-843).  The  basic 
inertial  navigation  system  is  described  by  a  nine  dimensional 
state  vector  X  consisting  of  position,  velocity,  and  air¬ 
craft  attitude  angles  with  respect  to  the  navigation  frame. 

An  east-north-up  orientation  of  the  navigation  frame  is  used 

in  the  IGI  simulator. 

* 
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X 
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\ 

Aircraft  Longitude 

L 

Aircraft  Latitude 

h 

Aircraft  Altitude 

Ve 

Aircraft  East  Velocity 

VN 

= 

Aircraft  North  Velocity  ( 

vz 

Aircraft  Up  Velocity 

€e 

Aircraft  Angular  Orientation  about  East  axis 

€N 

Aircraft  Angular  Orientation  about  North  axis 

_€z_ 

Aircraft  Angular  Orientation  about  Up  axis 

Evolution  of  the  states  is  described  by  the  following  system 
of  nonlinear  differential  equations 

X  =  F(X,  f,  t)  +  u(t)  (33) 

where  f  is  the  specific  force  and  u(t)  are  the  continuous  or 
discontinuous  rates  of  correction  in  an  aided-inertial  imple 
mentation.  For  example,  vertical  channel  stabilization  is 
accomplished  by  altitude  and  vertical  velocity  rates  from 
the  aircraft  altimeter. 

In  lieu  of  integrating  the  nonlinear  equations  in  (33) 
for  each  IGI  simulation  run,  the  approach  is  to  linearize 
equation  (33)  about  the  specified  (and  precomputed)  non¬ 
linear  flight  trajectory.  The  errors  are  propagated  as  per¬ 
turbations  to  the  reference  trajectory  and  added  to  the 
whole-valued  reference  trajectory  quantities  in  succeeding 
runs  to  obtain  a  first  order  approximation  to  the  actual 


inertial  navigation  system  outputs.  The  true  nonlinear 


flight  trajectory  consists  of  the  vector  Xt  (output  from 
PROFGEN  is  expressed  in  a  level  wander  angle  coordinate 
frame,  as  opposed  to  the  East-North-Up  geographic  frame) 
(See  Figure  11 ) . 


True  Longitude 

Lt 

True  Latitude 

ht 

True  Altitude 

Vet 

True  East  Velocity 

VNt 

True  North  Velocity 

vzt 

True  Up  Velocity 

fet 

= 

True  East  Specific  Force 

fNt 

True  North  Specific  Force 

Fzt 

True  Up  Specific  Force 

at 

True  Wander  Angle 

*t 

True  Aircraft  Roll  Angle 

9t 

True  Aircraft  Pitch  Angle 

It_ 

True  Aircraft  Yaw  Angle 

The  linearized  errors  SX^.  are  computed  from  equations  (33) 
and  (34)  using  the  following  set  of  differential  equations 


Figure  ||  Definition  of  PROFGEN  Geographic  <XG,YG,ZG) 
and  Wander  Angle  (X. ,Y*,ZJ  Coordinate  Frame  (17:7) 


where 


At(t)  =  (6F/6X) 


(36) 


are  the  partial  derivatives  of  F  (from  equation  (33))  eval¬ 
uated  on  Xfc.  The  w(t)  and  W*(t)  represent  stochastic  error 
forcing  functions  which  perturb  the  linear  system  dynamics 
(e.g.  gravity  anomalies,  gyro  biases,  and  accelerometer 
errors) . 

The  IGI  direct  simulation  is  currently  configured  for 
an  altimeter-aided  LN-15  inertial  navigation  system  error 
model.  The  functions  w(t),  w1(t),  and  u(t)  for  this  system 
are  derived  in  Reference  18,  The  w(t)  can  be  written  as 
follows. 


w(t )  = 


w 

■=R 

w„. 

= 

— € 

0 

Cp5fP 


,N 


+  6g 
P 


Position  Error  Components 
Velocity  Error  Components 
Attitude  Error  Components 


,N  . 


where  Cp  is  the  small  angle  transformation  matrix  from  plat- 


P  . 


form  to  navigation  coordinates  (East-North-Up);  6_f  is  the 


specific  force  measurement  error  in  platform  coordinates; 

1N  .  ... 

6g  is  the  gravity  model  error  in  navigation  coordinates; 


and  6w^  is  the  platform  angular  velocity  error  in  inertial 


coordinates .  These  INS  error  models  contain  a  total  of  48 


states  in  the  referenceu  report  (17»16).  In  addition  to  the 
48  states  from  the  INS  error  model,  seven  states  have  been 


added  to  simulate  the  GPS  UE  crystal  clock  errors,  two  ad¬ 
ditional  barometric  altimeter  errors  have  been  added,  and 
four  external  pseudo-range  error  states  have  been  added  (to 
account  for  the  small  inaccuracies  in  the  GPS  Control  and 
Space  Segments)  for  a  total  truth  state  vector  dimension  of 
61.  A  list  of  the  total  truth  states  is  provided  in  Table 
II  (17i6-50).  The  time  propagation  of  the  filter  states  is 
performed  with  the  upper  left  16  x  16  portion  of  the  Afc(t) 
matrix.  This  upper  left  matrix  is  referred  to  as  A^t)  and 
is  generated  using  nominal  position,  velocity  and  specific 
force.  Since  specific  force  is  not  directly  measurable  in 
an  inertial  system  an  approximation  based  on  velocity  dif¬ 
ferencing  is  used  (16j844).  The  generation  of  Af(t)  as  just 
described,  serves  to  isolate  the  user’s  inertial  system  state 
estimates  from  the  truth  state  values,  lending  added  realism 
to  the  simulation. 

IGI  Tra iectory  and  Navigation  Error  Formulation  (16: 
843-845) .  A  number  of  trajectories  and  navigation  errors 
are  employed  in  the  IGI  simulation  to  create  and  maintain 
the  distinction  between  the  user's  inertial  system  estimates 
and  the  truth  state  variables; 

(1)  True  Trajectory,  Xt 

(2)  True  Navigation  Errors,  5Xfc 

(3)  Nominal  Trajectory,  X^ 

(4)  Filter  Estimated  Navigation  Errors,  *Xf 

(5)  Filter  Estimated  Trajectory,  X- 


The  true  trajectory,  Xfc  is  generated  a  priori  by  the 
program  PROFGEN  (Ref.  24).  PROFGEN  processes  realistic 
flight  profiles  (path  accelerations,  altitude  rates,  head¬ 
ings,  etc.)  to  create  a  data  file  of  the  true  aircraft  state 
Xt  as  a  function  of  time.  The  processing  of  the  realistic 
flight  profile  input  is  accomplished  by  integrating  the  non¬ 
linear  equations  of  motion  over  a  WGS-72  (World  Geodectic 
Survey-1972)  reference  ellipsoid  with  a  Kutta-Merson  fifth- 
order  numerical  integration  routine  with  step-size  control 
(Ref.  24). 

The  true  nonlinear  reference  data  from  PROFGEN  is  used 
to  "drive”  the  linear  error  model  in  equation  (35).  Note 
that  the  true  specific  force  f  is  employed  directly  in  Afc(t) 
to  simulate  the  true  acceleration  "felt"  by  the  accelero¬ 
meters.  The  stochastic  disturbances  in  w^t)  are  generated 
by  means  of  a  Gaussian  random  number  generator  with  vari¬ 
ances  as  specified  in  Reference  18.  Since  the  steps ize  in 
equation  (35)  is  small,  typically  0.1  seconds,  integration 
of  the  equations  for  the  state  transition  matrix  t^ 

is  carried  out  with  the  linear  approximation. 


(38) 


and  the  resultant  true  error  is  given  by 


t . )6XJ  +  W1L 


6X  *+1  =  I_(t.^, 


(39) 


where  E[wlL]  =  0  ,  E[w1:LW1:iT]  =  flf(ti)6ij  and  W11  are  the 
integrated  contributions  of  the  noise  W1  in  analytical  form. 
The  true  initial  condition  errors  SX0.^  and  the  noise  covar¬ 
iance  are  specified  in  the  input  initialization  and  con¬ 
figuration  file. 

The  nominal  trajectory  which  represents  the  user’s 
INS-indicated  position  and  velocity,  is  created  by  the  dif¬ 
ference  between  the  true  navigation  error  6Xt  and  the  true 
reference  trajectory  X*. . 


*n  *  St  -  S*t  <40) 

Without  external  measurement  data,  or  between  measurement  up- 
dates,  Xjj  represents  the  user's  best  estimate  of  his  naviga- 
tion  parameters  based  solely  on  the  user’s  internal  naviga¬ 
tion  aids;  namely  the  IMU,  altimeter,  and  clock. 

The  nominal  nonlinear  reference  trajectory  X^  is  used 
to  "drive"  the  user's  16  state  filter  error  model.  The  in¬ 
puts  required  for  the  fundamental  matrix  A^(t)  include  nomi¬ 
nal  position,  velocity,  and  specific  force  f^.  Since  the 
specific  force  f^  is  not  directly  measurable  in  an  inertial 
system,  velocity  pulses  (or  differences)  are  utilized  to 
derive  an  approximation  as  follows 

fj^t)  =  [VN(t)  -  V^t  -  At)]/At  (41) 

where  At  is  the  velocity  pulse  integration  stepsize.  This 
duplicates  the  procedure  employed  in  a  real  inertial  system. 
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Therefore,  the  effects  of  integration  errors  and  time  lags 
(due  to  the  stepsize)  in  obtaining  correct  specific  force 
data  from  the  computer  will  be  reflected  in  the  simulation. 
Since  At  is  small,  the  following  approximation  is  made  for 
the  user's  estimate  of  the  state  transition  matrix. 


ti)  = 


i  + 


-  V 


(42) 


and  the  user's  estimated  navigation  error  is 

s4+1  =  ^f(ti+i*  v  6^f +  ^f1  (43) 

At  the  initial  time,  5X°  =  0  ,  since  the  user  has  no  better 
information  than  has  initial  nonlinear  state  estimate.  When 
satellite  measurements  are  available,  the  error  estimate 
6Xj.  is  updated  (measurement  update  will  be  discussed  further 
later  in  this  section) . 

The  user's  best  estimate  of  his  nonlinear  trajectory 
is  obtained  by  adding  the  error  estimate  6Xf  from  the  filter 
to  the  nominal  reference  trajectory, 

Xf  -  £n  +  (44) 

which  is  referred  to  as  the  estimated  trajectory.  The  true 
navigation  error  can  then  be  computed  by  talcing  the  differ¬ 
ence  between  the  estimated  and  true  reference  trajectories 


as  follows, 


or  using  equations  (40)  and  (44), 

=  SX^  -  5Xfc  (46) 

Satellite  Measurement  Generation  and  Incorporation  (17: 
42-50) .  The  IGI  simulator  assumes  that  the  18  GPS  satellites 
are  point  masses  moving  in  circular  orbits  of  12  hour  period 
above  a  spherical  earth.  Two-body  orbit  equations  and  angu¬ 
lar  transformations  are  used  to  compute  the  position  Xg  and 
velocity  of  the  selected  satellites.  The  position  and 
velocity  components  in  earth  centered  earth  fixed  (ECEF)  co¬ 
ordinates  at  a  given  time  are  computed  according  to  the  fol¬ 
lowing  equations  (See  Ref.  29:51-143). 

X  .  =  R„  sin  W  sin  i 
si  s 

X  -,  =  -  R  (sin  W  cos  i  cos  0  +  cos  W  sin  Q)  (47) 

s2  s 

X  -  =  -  R  (sin  W  cos  i  sin  Q  -  cos  W  cos  Q) 

Sj  S 

V  ,  =  V  cos  W  sin  i 

si  s 

v  =  -  V  (cos  W  cos  i  cos  fi  -  sin  W  sin  D)  +  X  - 

SZ  S  6  S  j 

V  ,  =  -  V  (cos  W  cos  i  sin  0  +  sin  W  cos  Q)  +  Q  X  7 

So  S  0  s z 

where 

R  =  radius  from  earth  center  to  the  satellite  orbit 
s 

V  =  orbital  velocity 


i  =  orbital  inclination 

W  =  argument  of  latitude  (See  Figure  12) 

Q  =  nodal  longitude  (See  Figure  12) 

0  =  earth  rotation  rate 

e 

The  ECEF  coordinate  frame  is  defined  by  an  axis  that  is  along 
the  earth’s  axis  of  rotation  (through  the  north  pole),  an 

.  *  .  ,  O 

axis  in  the  equatorial  plane  at  a  longitude  of  -90  from 
Greenwich,  and  the  third  axis  in  the  equational  plane  at  a 
longitude  equal  to  the  Greenwich  meridian.  The  simulated 
truth  state  satellite  pseudo-range  measurement  is  defined 
by  the  following 


PRm(t)  =  Rt(t)  +  6tu(t)  +  6t0(t)  +  Epr  (49) 


where 


PR  (t) 
m 

Rt(t) 


6te(t) 


Measured  pseudo-range 

Geometric  Range  from  the  user-to-satellite 
based  on  the  true  user’s  position 
True  user-clock  phase  error 

External  pseudo-range  error  due  to  satellite 
clock  errors,  satellite  ephemeaides  errors, 
and  atmospheric  delays 

Pseudo-range  additive  random  error  (due  to 
multipath  random  error,  receiver  noise,  etc.) 


t  =  Pseudo-range  measurement  time 


y1 


vernal  equinox 
direction 


satellite  *  position 


where  s 

Periapsis  Direction  =  the  direction  of  closest  orbital 

approach  to  the  earth 

h  =  the  specific  angular  momentum  =  r  X  V 

e  =  eccentricity  vector 

I,  J,  and  K  denote  the  Geocentric- 
Eguational  Coordinate  System 

vQ  =  True  Anomaly  at  epoch 

n  =  Longitude  of  the  ascending  node 
00  =  Argument  of  periapsis 

Figure  12.  Classical  Orbital  Element  Description  (29s 59). 
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The  clock  phase  errors  6t  and  external  pseudo-range 

error  6t  are  included  as  states  in  the  truth  model  (17: 35— 

36).  The  pseudo-range  additive  random  error  E  is  gener- 

ated  by  a  random  number  generator  (Gaussian  distributed  with 

zero  mean  and  unit  variance)  times  the  one-sigma  value  of  a 

pseudo-range  additive  random  error  (7  feet).  The  simulated 

filter  state  pseudo-range  measurement,  however,  does  not 

have  6t  and  6t  available  in  the  filter  states.  To  account 

for  the  added  uncertainty  in  the  simulated  filter  state 

pseudo-range  due  to  the  lack  of  knowledge  of  6t  and  6t  , 

the  one-sigma  value  of  pseudo-range  additive  random  error 

E  has  a  value  of  12  feet  (17?44). 
pr 

Measurement  information  is  processed  in  the  Kalman  fil¬ 
ter  by  using  the  measurement  residual.  The  measurement  re¬ 
sidual  is  calculated  by  differencing  the  pseudo-range  meas¬ 
urement  with  the  expected  value  for  the  pseudo-range  con¬ 
ditioned  on  the  previous  measurements  taken.  The  expected 
value  is  computed  based  on  the  INS -indicated  position  cor¬ 
rected  by  the  Kalman  filter  estimates  of  error  in  the  indi- 

A  A  A 

cated  position  (6Xf,  6Lf,  6hf)  and  the  given  satellite 
ephemeris.  The  equation  for  the  expected  value  of  the  pseudo 
range  measurement  is  (the  detailed  development  of  equation 
(50)  is  available  in  Reference  17,  page  45-48), 

PRf(t)  =  gpR(rN(t),  5Xf(t))  +  hp^SXj (t))E[(SXf (t) 

-  6X-(T))]  +  E  [V  ] 


(50a) 


f  = 


for  i  =  1  -•  m 
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=  o 


for  j  =  1  -*-Nf 
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j 
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sfj<+> 


■  Csf j(_)  - 


-  ^N^aN 


A  A 

6Xf (  +  )  =  63^ (-)  +  K  A  Zi 
Pf(+)  =  Sf(+)SfT(+) 


(53a) 

(53b) 

(53c) 

(53d) 

(53e) 

( 53f ) 

(53g) 

(53h) 

(54) 


where 


Pf(*)  =  filter  covariance  matrix  before  (-)  or  after 
(+)  measurement  incorporation 
Sf( •)  =  the  upper  triangular  square  root  matrix  of  the 

filter  covariance  (using  Cholesky  decomposition) 
Sf  j  =  the  j*'*1  column  of  Sf 
f j  =  the  j  element  of  f  (specific  force) 

AZ^  =  the  i  measurement  residual 


6 


=  the  measurement  gradient  vector  of  the  i 
measurement 

2  tli. 

a ^  =  variance  of  the  noise  in  the  i  measurement 

Nf  =  dimension  of  the  filter  state 
M  =  number  of  measurements 
K  =  Kalman  gain  vector 

Modif ications  to  the  IGI  Simulator 

The  modifications  to  the  IGI  Simulator  necessary  to 
test  the  cost  criterion  satellite  selection  algorithm  fall 
into  two  categories,  functions  supporting  the  cost  criterion 
and  functions  to  increase  the  realism  of  the  computer  simula¬ 
tion.  The  increased  realism  will  allow  the  results  of  this 
study  to  be  more  applicable  to  the  present  hardware  design. 
Also,  the  increased  realism  and  capability  of  the  modified 
IGI  simulator  will  make  it  an  effective  design/analysis  tool. 
The  major  functions  added  to  the  IGI  simulator  are:  the  cost 
criterion  (for  satellite  selection),  an  adaptive  null  steer¬ 
ing  antenna  algorithm,  the  air  and/or  ground  jammers,  and  an 
adaptive  bandwidth  control  of  the  GPS  satellite  signal  track¬ 
ing  loops.  A  brief  discussion  of  each  function  follows, 
with  references  to  the  appropriate  appendix  for  a  more  de¬ 
tailed  analysis  of  each  function. 

The  Cost  Criterion.  The  cost  criterion  function  added 
to  the  IGI  simulator  modifies  the  satellite  selection  algo¬ 
rithm  to  minimize  equation  (33)  instead  of  GDOP  for  each  pos¬ 
sible  satellite  set.  In  order  to  calculate  equation  (33), 
the  covariance  matrix  values  for  three  position  errors  and 


one  time  bias  error  are  needed  to  fill  the  a  priori  covari¬ 
ance  matrix  P^.  Also,  the  one-sigma  measurement  noise  values 
for  each  satellite  (in  the  set  being  examined)  are  needed  to 
create  the  measurement  noise  matrix  R.  The  information 
needed  for  the  matrix  is  readily  available  from  the  IGI 
simulator.  The  information  for  the  R  matrix,  however,  is  a 
function  of  satellite  signal  level  at  the  antenna,  the 
jamming  power  level  at  the  antenna,  the  type  of  antenna,  and 
the  type  of  satellite  tracking  loop,  none  of  which  are  avail¬ 
able  in  the  IGI  simulator,  but  are  added  in  the  modifications 
to  be  discussed.  For  analysis  purposes,  the  cost  criterion 
value  of  the  satellite  set  with  the  minimum  cost  criterion 
value  is  stored  for  output  with  the  printed  data  as  well  as 
the  GDOP  value  for  that  particular  satellite  set.  Then  by 
comparing  the  cost  and  GDOP  values  for  a  simulation  run 
using  the  cost  criterion  and  a  simulation  run  using  GDOP, 
conclusive  results  can  be  made  of  the  effectiveness  of  the 
alternate  versus  standard  satellite  selection  criterion  on 
mean  square  position  error. 

The  Null  Steering  Antenna .  The  objective  of  the  null 
steering  antenna  algorithm  is  to  maximize  the  GPS  satellite 
signals,  minimize  interference,  and  estimate  the  signal 
carrier-to-noise  density  (C/No)  values  for  each  in  view  sat¬ 
ellite.  To  accomplish  these  objectives  the  null  steering 
algorithm  uses  a  power  inversion  process  to  determine  the 
directions  of  the  strongest  interference  sources.  The  null 
steering  controller  thereby  forms  nulls  in  the  direction  of 


the  interference  sources  to  reduce  their  jamming  power.  The 
implemented  null  steering  an'  -?nna  model  uses  a  seven  element 
phased  array  antenna  arranged  in  a  circular  pattern  (Figure 
13).  A  seven  element  null  steering  antenna  is  capable  of 
forming  a  maximum  of  six  nulls  with  approximately  30db  watts 
of  jammer  suppression  (6:141).  The  satellite  signals  are 
protected  from  spurious  nulls  by  using  the  estimate  of  the 
satellites  direction  in  the  null  steering  controller. 

The  estimate  of  signal-to-noise  density  (C/No)  values 
for  each  satellite  is  an  inherent  quantity  in  the  null  steer¬ 
ing  algorithm  which  limits  the  unwanted  signals  and  maxi¬ 
mizes  the  satellite  signals. 

For  a  more  detailed  analysis  of  the  processes  implemented 
in  the  modified  IGI  simulator,  the  reader  is  directed  to 
Appendix  B. 

Jamming  Field.  The  jamming  field  parameters  are  input 
through  the  initialization  and  configuration  input  file.  The 
simulation  user  may  specify  the  number  of  jammers  (1-25), 
the  position  of  the  jammers  (in  degrees  latitude,  degrees 
longitude,  and  altitude  in  feet),  the  velocity  of  the  jammers 
(in  feet-per-second) ,  and  the  effective  isotropic  radiating 
power  (EIRP)  out  of  the  jammer.  The  above  format  was  imple¬ 
mented  to  allow  land,  sea,  and  airborne  jammers  that  are 
stationary  or  mobile.  The  jammers  are  omni-directional 
radiating  sources  that  transmit  the  full  EIRP  into  the  GPS 
LI  frequency.  The  GPS  user-to- jammer  line-of -sight  is 


Right  Wing  of 
the  Aircraft 


Figure  13.  Seven  Element  Null  Steering  Antenna  Geometry 


obstructed  by  blockage  from  the  earth  only  (which  is  again 
assumed  to  be  spherical)  (See  Appendix  B) . 

Adaptive  Bandwidth  Control  of  Satellite  Tracking .  As 
discussed  in  Chapter  III,  the  bandwidth  of  the  satellite 
tracking  loops  is  very  important  in  limiting  jamming  related 
navigation  errors.  The  C/No  estimates  from  the  null  steering 
algorithm  allow  the  implementation  of  adaptive  control  of 
the  bandwidth  of  the  satellite  tracking  loops.  The  band¬ 
width  of  the  code  loops  is  determined  by  four  separate  quad¬ 
ratic  functions  that  use  the  C/No  estimate.  The  software 
implementation  covers  both  coherent  and  non-coherent  code 
loop  modes,  as  well  as  loss-of-lock  thresholds.  The  carrier 
loop  was  not  implemented  due  to  the  very  small  amount  of 
error  associated  with  the  carrier  loop  estimates  (~  1  cm  of 
position  error)  (5:25-26).  A  more  detailed  discussion  of 
the  adaptive  bandwidth  control  of  satellite  tracking  is  ad¬ 
dressed  in  Appendix  B. 


Introduction 


This  chapter  will  present  the  results  from  the  modi¬ 
fied  IGI  computer  simulation  runs.  The  results  will  be  pre 
sented  in  a  pictorial  and  tabular  form  to  increase  readabil 
ity.  From  the  abundance  of  IGI  output  information  four 
parameters  will  be  highlighted! 

(1)  The  number  of  times  the  cost  criterion  chooses  a 
satellite  set  other  than  the  set  chosen  using  GDOP. 

(2)  The  user  position  and  clock  phase  errors  for  both 
GDOP  and  the  cost  criterion  for  the  times  when  an  alternate 
satellite  set  was  chosen  from  GDOP. 

(3)  The  measurement  error  estimate  for  each  satellite 
at  the  satellite  selection  time. 

(4)  The  noise  level  at  each  satellite  selection  time 
when  GDOP  and  the  cost  criterion  don't  agree  on  the  same 
satellite  set. 

These  four  parameters  will  be  the  basis  for  an  assessment 
of  the  performance  of  GDOP  versus  the  cost  criterion  in 
satellite  selection. 

The  Simulated  Mission 

The  mission  used  in  all  of  the  simulation  runs  is  an 
F-4  Close  Air  Support  (CAS)  flight  profile,  provided  by  the 
Reference  Systems  Branch,  AFWAL  Avionics  Laboratory,  Wright 
Patterson  AFB,  Ohio.  This  flight  profile  was  chosen  be¬ 
cause  of  the  high  aircraft  dynamics  associated  with  the 


tactical  mission  and  the  increased  realism  of  the  mission 
by  locating  a  jamming  field  around  the  target  area.  A  24 
minute  segment  of  the  total  F4  CAS  mission  was  used  in  the 
simulation  runs  in  order  to  cover  the  aircraft’s  maneuver’s 
in  the  target  area,  while  keeping  the  simulation  outputs  to 
a  manageable  level.  The  time  history  of  the  24  minute  F4 
CAS  mission  is  presented  in  Table  III.  The  target  is  as¬ 
sumed  to  be  located  at  the  position  where  the  F-4  is  at  its 
lowest  altitude  (time  =  760.0  seconds). 

The  jamming  field  patterns  for  each  simulation  run 
(simulation  runs  are  identified  as  missions  1-5)  were  modeled 
as  stylized  jamming  arrays  (as  in  Ref.  20)  to  easily  vary 
the  maximum  jamming  power,  to  move  the  maximum  jamming  power 
points,  and  to  increase  or  decrease  the  rate  of  change  of 
jamming  power.  The  location  and  transmitted  power  for  each 
jammer  in  each  mission  in  shown  in  Table  IV. 

Mission  Parameters 

Each  simulated  mission  has  a  duration  of  1400  seconds. 

At  time  =  0.0  seconds  the  truth  state  terms  and  the  filter 
covariance  terms  of  east,  north,  and  up  user  position  error 
and  the  clock  phase  error  are  set  to  a  value  of  50  feet  (that 
is  a  50  ft  one-sigma  for  the  covariance  terms).  This  is  con¬ 
sistent  with  an  integrated  GPS/INS  that  is  within  its  steady 
state  performance  accuracy,  attained  enroute  to  the  target 
area.  The  GPS  measurements  are  incorporated  every  five 
seconds,  but  the  truth  and  filter  state's  and  the  filter  co- 
variance  data  is  printed  out  once  every  100  seconds.  While 
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Jamming  Field  Patterns 
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All  jammers  have  an  elevation  of  10  feet  and  are  station 
ary  with  respect  to  the  earth's  surface.. 
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in  the  jamming  field,  satellite  selection  using  both  the  GDOP 
and  cost  criterion  is  executed  every  100  seconds.  There  is 
no  measurement  incorporation  or  satellite  selection  at  the 
final  time  (1400  seconds). 

Mission  #  1_ 

The  satellite  selection  in  Mission  #  1  between  GDOP  and 
the  cost  criterion  was  in  agreement  until  700  seconds  into 
the  flight.  For  next  6  satellite  selection  times,  a  set 
other  than  the  best  GDOP  set  was  chosen  by  the  cost  criterion 
Table  Va  shows  the  measurement  noise  estimates  for  the  satel¬ 
lites  chosen  by  the  cost  criterion  and  the  associated  noise 
density  at  the  6  satellite  selection  times  of  interest.  The 
noise  density  (No)  value  is  N^,  +  Nj  from  Equation  (74)  in 
Appendix  B.  The  Nj  term  is  the  difference  of  the  total 
jammer  power  at  the  antenna  output  ( J  )  and  the  p-code  proc- 

a 

essing  gain  (6)  when  both  J_  and  G  are  expressed  in  decibels. 
The  noise  density  (No)  is  related  to  another  commonly  used 
term  J/S.  The  term  J.S  is  defined  as  the  jammer -to-signal 
power  ratio  (typically  expressed  in  dBs).  A  noise  density. 

No  =  190  db  relates  to  J/S  values  of  45-58.  J/S  values  are 
typically  associated  with  satellite  signal  and  jammer  power 
estimates  with  respect  to  the  antenna  input  (i.e.,  the  anti¬ 
jam  performance  gained  from  a  null  steering  antenna,  adaptive 
bandwidth  control,  and  the  p-code  processing  gain  is  lost), 
therefore,  the  C/No  or  No  terms  provide  much  more  insight 
into  the  UE  sets'  performance  with  respect  to  jamming  power. 


Table  Vb  presents  the  same  data  as  Table  Va  for  the  GDOP 
selected  satellite  set.  Note  that  the  cost  criterion  will 
almost  always  choose  satellites  that  have  measurement  errors 
that  are  less  than  the  measurement  errors  for  the  GDOP  se¬ 
lected  satellites.  The  performance  of  the  integrated  GPS/INS 
with  the  cost  criterion  versus  the  GDOP  satellite  selection 
algorithm  is  presented  in  Figures  14-17.  broken  into  the 
three  components  of  filter  estimation  (filter-truth  states) 
position  error  and  filter  estimation  clock  phase  error  (note 
that  plot  a  in  each  figure  is  the  IGI  using  the  cost  crite¬ 
rion,  and  plot  b  in  each  figure  is  the  IGI  using  GDOP).  In¬ 
cluded  in  the  plots  of  Figures  14-17  are  the  ±  1  sigma 
values  of  the  related  filter  covariance  terms.  The  figures 
show  the  cost  criterion  exhibiting  smaller  errors  in  all 
three  components  of  position  error  than  the  GDOP  plots,  with 
the  clock  phase  errors  being  roughly  equal  during  the  700- 
1200  second  interval  when  an  alternate  satellite  set  is 
being  tracked.  The  maximum  noise  density  is  -186.93dB  and 
occurs  at  approximately  1053  seconds.  The  composite  statis¬ 
tics  for  the  filter  estimation  radial  position  errors  aver¬ 
aged  over  the  700-1200  second  interval  arej 

Cost  Criterion  versus  GDOP 

Mean  Std.  Dev.  Mean  Std.  Dev. 

241.77  262.38  627.80  855.04 

Including  time-bias  range  errors i 

Mean  Std.  Dev.  Mean  Std.  Dev. 

294.81  262.38  627.80  855.04 

where  all  units  are  in  feet. 
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Figure  14a.  Mission  #  li  Cost  Criterion;  East  Position  Error 


GDOPj  East  Position  Error 


Mission 


Figure  16b.  Mission  #  li  GDOP;  Altitude  Error 


Figure  17a.  Mission  #  Is  Cost  Criterion;  Clock  Phase  Error 


Mission  #  2 

The  satellite  selection  in  Mission  #  2  varied  from 
700-1100  seconds  and  1200-1300  seconds  in  the  flight.  Table 
Via  and  b  shows  measurement  error  estimates  for  the  cost 
criterion  and  GDOP  satellite  sets  respectively.  Note  that 
the  noise  density  peaks  much  earlier  in  Mission  #  2  than  in 
Mission  #  1  (max  noise  density  for  Mission  #  2  is  -179.24  dB 
and  occurs  at  747  seconds  in  the  mission),  then  remains  rel¬ 
atively  low  for  the  duration  of  the  mission.  As  was  the  case 
in  Mission  #  1,  the  cost  criterion  has  chosen  a  satellite 
set  with  a  lower  composite  measurement  error  than  the  GDOP 
satellite  set. 

Figures  18-21  present  the  time  history  performance  of 
position  and  clock  phase  errors.  The  figures  show  that  the 
cost  criterion  exhibits  smaller  position  and  clock  phase 
errors  than  the  corresponding  GDOP  errors,  but  that  the  dif¬ 
ference  is  much  less  than  appears  in  Mission  #  1.  The  com¬ 
posite  statistics  for  the  filter  estimation  radial  position 
errors  averaged  over  the  700-1300  second  interval  (excluding 
1100  seconds)  are: 

Cost  Criterion  versus  GDOP 

Mean  Std.  Dev.  Mean  Std.  Dev. 

73.81  34.35  94.20  45.13 

Including  time-bias  range  errors: 

Mean  Std.  Dev.  Mean  Std.  Dev. 

91.12  34.46  116.13  38.69 

where  all  units  are  in  feet.  Again,  the  cost  criterion  sig¬ 
nificantly  out-performs  the  GDOP  criterion. 
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Figure  18a.  Mission  #  2i  Cost  Criterion;  East  Position  Error 


Figure  20b.  Mission  #  2i  GDOP;  Altitude  Error 
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Figure  21a.  Mission  #  2:  Cost  Criterion;  Clock  Phase  Error 


Mission  #  3 

The  satellite  selection  using  the  cost  criterion  varied 
from  the  GDOP  choice  at  700  seconds  and  for  the  following  6 
satellite  selection  times.  Table  Vila  and  b  shows  the  mea¬ 
surement  errors  for  the  cost  criterion  and  GDOP  satellite 
sets  respectively. 

The  cost  criterion  chose  a  satellite  set  with  a  lower 
composite  measurement  error  than  the  GDOP  satellite  set  for 
every  selection  time  except  at  1300  seconds.  The  maximum 
noise  density  is  -194.50dB  at  747  seconds.  Note  the  slow 
changes  in  the  noise  density  from  one  satellite  selection 
time  to  the  next.  Also,  the  measurement  error  estimates  for 
the  four  satellite  sets  are  roughly  equivalent  (i.e.,  Rmeas 

<!>  ~  Rmeas<2>  »  RmeaS<3)  ~  Rmeas‘4>-  This  effect  TOUld 
make  the  R  matrix  approximately  a  scalar  multiple  of  the 

identity  matrix.  As  stated  in  Chapter  III  for  noise  weighted 
GDOP,  if  the  R  matrix  is  the  identity  matrix  or  a  scalar 
multiple  of  the  identity  matrix  then  the  satellite  selection 
criterion  is  the  same  as  GDOP  values  but  increased  by  the 
scalar  multiple.  Applying  this  effect  to  the  cost  criterion 
means  that  there  is  a  trade-off  between  the  GDOP  effect  and 
the  filter  covariance's  need  for  information  in  a  particular 
direction.  Since  the  satellite  set  chosen  by  the  cost  in¬ 
terior  differs  from  the  GDOP  selection,  the  cause  can  be  at¬ 
tributed  to  the  filter  covariance  terms  outweighing  the  GDOP 
terms . 

Figures  22-25  present  the  time  history  performance  of 
position  and  clock  phase  errors  for  this  mission.  The 
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22a.  Mission  #  3j  Cost  Criterion;  East  Position  Error 


Figure  23a.  Mission  #  3:  Cost  Criterion;  North  Position  Error 


Figure  23b.  Mission  #  3:  GDOPj  North  Position  Error 


Figure  24a.  Mission  #  3i  Cost  Criterion;  Altitude  Error 


Figure  24b.  Mission  #  3t  GDOP;  Altitude  Error 


Figure  25a.  Mission  #  3i  Cost  Criterion;  Clock  Phase  Error 


Figure  25b.  Mission  #  3i  GDOP;  Clock  Phase  Error 


figures  show  a  small  improvement  in  performance  for  the  cost 
criterion  over  the  GDOP  case.  The  composite  statistics  for 
the  filter  estimation  radial  position  errors  averaged  over 
the  700-1300  second  interval  are: 

Cost  Criterion  versus  GDOP 

Mean  Std.  Dev.  Mean  Std.  Dev. 

74.88  42.23  85.03  37.95 

Including  the  time-bias  range  errors: 

Mean  S td .  Dev .  Mean  Std.  Dev. 

89.49  40.33  107.07  37.96 

where  all  units  are  in  feet.  The  cost  criterion  shows  a 
small  performance  improvement  over  GDOP,  but  the  larger 
standard  deviation  of  the  cost  criterion  would  reduce  the  sig¬ 
nificance  of  the  difference  in  mean  values. 

Mission  #  4 

As  in  Mission  #  3,  the  satellite  set  chosen  using  the 
cost  criterion  varied  from  the  GDOP  choice  from  700  seconds 
to  1300  seconds  in  the  mission.  The  maximum  noise  density 
is  -172.40db  at  740  seconds  into  the  mission.  Table  Villa 
and  b  shows  the  measurement  errors  for  the  cost  criterion 
and  GDOP  satellite  sets  respectively.  The  cost  criterion 
chose  a  satellite  set  with  larger  composite  measurement  error 
than  the  GDOP  satellite  set  at  the  800  and  900  second  selec¬ 
tion  times .  This  choice  may  be  a  result  of  the  rapid  change 
in  the  noise  density  and  encountering  the  maximum  noise  den¬ 
sity  unexpectedly  (i.e.,  without  some  significant  build-up 
in  noise  density  at  700  seconds). 
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Figures  26-29  present  the  time  history  performance  of 
the  position  and  clock  phase  errors  for  this  mission.  The 
figures  show  that  the  performance  is  roughly  comparable  be¬ 
tween  the  cost  criterion  and  GDOP.  The  cost  criterion  fil¬ 
ter  covariance  reaction  to  the  increased  noise  density  from 
700-960  seconds  probably  accounts  for  the  roughly  equivalent 
performance  errors  with  respect  to  GDOP,  even  though  a 
noisier  measurement  set  is  chosen  at  800  and  900  seconds. 

The  filter  covariance  reaction  to  increasing  or  decreasing 
noise  density  can  be  viewed  as  the  way  the  cost  criterion 
keeps  track  of  the  time  history  of  the  noise  density.  For 
instance,  small  a  priori  filter  covariance  values  with  large 
measurement  errors  would  indicate  entry  into  the  jamming 
field.  When  both  a  priori  filter  covariance  values  and 
measurement  errors  are  large  the  user  is  operating  within 
the  jamming  field.  When  the  a  priori  filter  covariances 
values  are  large  and  the  measurement  errors  are  small  the 
user  is  leaving  the  jamming  field.  How  effectively  the  cost 
criterion  uses  the  jamming  information  seems  to  be  dependent 
on  the  rate  of  change  of  the  noise  density.  The  slower  the 
variation  in  noise  density  the  more  significant  the  improve¬ 
ment  in  performance  for  the  cost  criterion  with  respect  to 
GDOP. 

The  composite  statistics  for  the  filter  estimation 
position  errors  for  the  700-1300  second  interval  are: 

Cost  Criterion  versus  GDOP 


Mean 


Std.  Dev 


Mean 


Std.  Dev. 


Figure  26a.  Mission  #  4:  Cost  Criterion;  East  Position  Error 


Figure  27a.  Mission  #  4:  Cost  Criterion;  North  Position  Error 


Figure  27b.  Mission  #  4s  GDOPj  North  Position  Error 


Figure  28a.  Mission  #  4:  Cost  Criterion*  Altitude  Error 


GDOP ;  Altitude  Error 


29a.  Mission  #  4s  Cost  Criterion;  Clock  Phase  Error 
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Including  the  time-bias  range  errors: 

Mean  Std.  Dev.  Mean  Std.  Dev. 

164.38  86.74  146.05  98.76 

where  all  units  are  in  feet.  In  Missions  1-3  the  time-bias 
range  errors  were  roughly  equal  for  both  the  cost  criterion 
and  GDOP,  or  they  did  not  significantly  affect  the  performance 
comparison.  In  this  mission  the  time-bias  range  errors 
change  the  performance  assessment  significantly.  The  large 
time  bias  range  errors  for  the  cost  criterion  can  be  ex- 
plained  by  the  (HR  H)  term  in  the  cost  criterion.  When  a 
satellite  set  with  poor  time-bias  related  geometry  has  a 
composite  measurement  error  much  smaller  than  a  more  favor¬ 
able  time-bias  geometry,  the  first  set  will  be  selected. 

Adding  the  filter  covariance  terms  into  the  selection  proc¬ 
ess  further  increases  the  desire  to  select  the  poor  time-bias 
related  geometry.  This  last  effect  is  due  to  the  fact  that 
under  the  influence  of  increased  noise  density,  the  filter 
covariance  position  terms  reflect  much  larger  errors  than 
does  the  filter  covariance  clock  phase  error.  The  direct  re¬ 
sult  in  the  cost  criterion  is  a  need  for  information  for 
position  errors  that  significantly  outweigh  the  clock  phase 
errors . 

Mission  #  5 

The  satellite  set  chosen  using  the  cost  criterion 
varied  from  the  GDOP  choice  during  700-1000  seconds  and 
1100-1300  seconds  in  the  mission.  The  maximum  noise  density 
was  -I63.60db  at  780  seconds  in  the  flight.  The  satellite 
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set's  measurement  errors  are  presented  in  Table  IXa  and  b. 
Note  that  at  80  seconds  the  GPS  signal  is  jammed-of f-the 
air  and  at  900  seconds  only  one  satellite  is  being  tracked 
for  either  the  cost  criterion  or  GDOP  case.  At  700,1100, 
1200,  and  1300  seconds  the  cost  criterion  has  selected  a 
satellite  set  with  lower  composite  measurement  error  than 
the  GDOP  selected  set.  The  encounter  with  the  maximum  noise 
density  results  in  over  a  30db  increase  in  noise  density  from 
700-800  seconds.  Figures  30-33  plot  the  time  history  per¬ 
formance  of  the  filter  estimation  position  and  clock  phase 
errors.  In  all  cases  the  ^DOP  performance  is  better  than 
the  cost  criterion.  The  composite  statistic  of  the  700-1300 
second  time  interval  ares 

Cost  Criterion  versus  GDOP 


Mean 

Std.  Dev. 

Mean 

Std.  Dev. 

125.59 

114.50 

87.84 

37.92 

Including  the 

time-bias  range 

errors : 

Mean 

Std.  Dev. 

Mean 

Std.  Dev. 

142.21 

110.01 

117.50 

39.34 

where  all  units  are  in  feet.  The  performance  degradation  of 
the  cost  criterion  with  respect  to  GDOP  is  even  more  pro¬ 
nounced  in  this  mission  than  in  Mission  #  4.  The  rapid  vari¬ 
ation  in  the  noise  density  is  also  more  pronounced  than  in 
Mission  #  4.  The  time-bias  range  errors  are  not  significant 
in  performance  assessment  in  this  mission  as  they  were  in 
Mission  #  4. 
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Figure  31a.  Mission  #  5:  Cost  Criterion;  North  Position  Error 


Figure  31b.  Mission  #  5j  GDOP;  North  Position  Error 
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32a 


GDOP;  Altitude  Error 


33a.  Mission  #  5s  Cost  Criterion;  Clock  Phase  Error 
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Figure  33b.  Mission  #  5s  GDOP;  Clock  Phase  Error 


Summary 

This  chapter  has  reviewed  the  data  from  5  simulated 
close-air-support  missions  with  an  integrated  GPS/INS  navi¬ 
gator  in  a  jamming  environment.  Each  mission  was  run  twice, 
once  using  the  cost  criterion  and  once  using  the  GDOP  cri¬ 
terion  for  satellite  selection.  Several  parameters  are  of 
particular  interest  in  the  assessment  of  the  performance  im¬ 
provement  of  the  cost  criterion  with  respect  to  the  GDOP 
method  for  satellite  selection.  Table  X  presents  a  review 
of  the  composite  statistics  for  all  five  missions.  It  is 
clear  from  Table  Xb  and  c  that  on  Missions  1  and  2  the  cost 
criterion  yields  a  performance  improvement  with  respect  to 
GDOP;  on  Mission  3  performance  is  equivalent  to  GDOP;  and 
Missions  4  and  5  indicate  performance  degradation  of  the  cost 
criterion  with  respect  to  GDOP.  The  performance  degrada¬ 
tion  of  the  cost  criterion  in  Missions  4  and  5  can  be  ex¬ 
plained  by  the  rapid  variation  in  noise  density  and  the  larger 
measurement  errors  of  the  selected  satellite  set  from  800  to 
1000  seconds.  This  would  suggest  that  the  cost  criterion 
would  have  to  be  tested  against  realistic  jamming  scenarios 
to  determine  its  value  for  a  particular  mission.  Note  that 
in  all  5  missions  for  the  first  600  seconds  only  5  satellites 
are  in-view,  but  during  the  next  700  seconds  of  the  mission 
(700-1400  seconds)  6  satellites  are  in-view.  This  accounts 
for  the  fact  that  all  five  missions  have  alternate  satellite 
sets  chosen  at  700  seconds  and  for  the  next  six  selection 


times  as  well  (with  the  exception  of  Missions  1,  2,  and  5 


TABLE  Xa 

Composite  Statistics 

No.  of  Alterna-  Max.  Noise  Time  of  Max. 


Mission  # _ tive  Sets _ Density _ Noise 


1. 

6 

out  of  14 

-186.93 

db 

watts 

1053  seconds 

2. 

6 

out  of  14 

-179.24 

db 

watts 

747  seconds 

3. 

7 

out  of  14 

-194.50 

db 

watts 

747  seconds 

4. 

7 

out  of  14 

-172.40 

db 

watts 

740  seconds 

5. 

6 

out  of  14 

-163.60 

db 

watts 

780  seconds 

0-700 

sec  =  5 

satellites 

700-1400  seconds  = 

6  satellites 

TABLE  Xb 

Performance  Comparison  of  User 

Position  Error 

(Without  Time-Bias  Range 

Errors) 

Cost 

Criterion 

versus 

GDOP 

Mean 

Std.  Dev 

• 

Mean 

Std.  Dev. 

1. 

241.77 

214.65 

604.09 

848.80 

2. 

73.81 

34.35 

94.20 

45.13 

3. 

74.88 

42.23 

85.03 

37.95 

4. 

120.31 

81.14 

123.64 

109.12 

5. 

125.59 

114.50 

87.83 

37.92 

1  all 

units  are  in  feet 

TABLE  Xc 


Performance  Comparison 
(With  Time-Bias 

of  User 
Range 

Position 

Errors) 

Error 

Cost 

Criterion 

versus 

GDOP 

Mean 

Std.  Dev. 

Mean 

Std.  Dev. 

1. 

294.81 

262.38 

627.80 

855.04 

2. 

91.12 

34.46 

116.13 

38.69 

3. 

89.49 

40.33 

107.07 

37.96 

4. 

164.38 

86.74 

146.05 

98.76 

5. 

142.21 

110.01 

117.50 

39.34 

/  all  units  are  in  feet 


that  chose  an  alternate  satellite  set  five  out  of  the  next 
six  selection  times).  The  limited  satellite  visibility  was 
a  concern  that  was  discussed  in  Chapter  Ill’s  analysis  of 
satellite  visibility  for  the  3  satellites  x  6  orbits  con¬ 
stellation.  Increasing  the  number  of  in-view  satellites 
(by  either  scheduling  missions  during  periods  of  high  visi¬ 
bility  or  by  increasing  the  number  of  active  satellites  in 
the  GPS  constellation)  increases  the  number  of  choices  avail¬ 
able  to  the  cost  criterion  and  the  performance  improvement 
over  the  GDOP  satellite  selection  criterion. 


VI  Conclusions/Recommendations 

Conclusions 

This  report  has  tested  and  analyzed  an  alternate  sat¬ 
ellite  selection  criterion  for  an  integrated  GPS/INS  naviga¬ 
tor  to  document  performance  improvements  with  respect  to  the 
standard  GDOP  criterion.  Testing  was  accomplished  using  a 
version  of  the  IGI  computer  simulation,  modified  specifically 
to  test  satellite  selection  criterion.  The  modifications  to 
the  simulation  included  the  addition  of  a  null  steering 
antenna  algorithm,  adaptive  bandwidth  control  of  the  satel¬ 
lite  tracking  loops,  the  parameters  that  describe  the  jamming 
field,  and  the  implementation  of  the  new  satellite  selection 
criterion  (known  as  the  cost  criterion) .  The  computer  sim¬ 
ulated  missions  were  created  by  running  a  single  F-4  Close- 
Air-Support  (CAS)  flight  profile  against  five  different  jam¬ 
ming  field  patterns  (where  the  power  of  the  jammer's  and/or 
the  dimensions  of  the  jamming  field  were  varied  from  mission 
to  mission).  The  results  from  the  five  missions  indicated 
that  the  performance  of  the  cost  criterion  was  dependent  on 
the  rate  of  change  in  noise  density  at  the  GPS  receiver  and 
the  number  of  satellites  available. 

For  a  slowly  varying  noise  density,  the  cost  criterion 
showed  performance  improvement  of  10-350  feet  in  mean  radial 
position  error  (with  respect  to  GDOP  performance) .  For 
rapidly  varying  noise  density  the  cost  criterion  showed  a 
performance  degradation  of  10-30  feet  in  mean  radial  position 
error  (with  respect  to  GDOP  performance) .  The  cost  criterion 


was  tested  in  Reference  11  (Brogan  Study)  using  a  static 
measurement  analysis  to  assess  the  performance  of  GPS  UE. 

The  static  versus  dynamic  measurement  analysis  looks  at  sat¬ 
ellite  measurement  selection  for  one  point  in  time  as  opposed 
to  satellite  measurement  selection  as  a  function  of  the  time 
history  (and  the  associated  variations  in  the  noise  density) 
of  the  mission.  The  performance  degradation  observed  in  this 
study  is  in  sharp  contrast  to  performance  improvements  sug¬ 
gested  by  the  static  measurement  analysis  of  the  cost  cri¬ 
terion  in  Reference  11.  The  reason  for  the  contradiction  in 
performance  between  the  static  and  dynamic  analysis  is  that 
the  static  measurement  case  made  no  attempt  to  account  for 
rapid  variations  in  noise  density.  During  rapid  variations 
in  the  noise  density,  the  cost  criterion  could  make  worse 
decisions  than  GDOP  (i.e.,  choose  a  satellite  set  resulting 
in  larger  position  errors  than  the  satellite  set  chosen  by 
GDOP)  due  to  satellite  selection  based  on  measurement  noise 
estimates  that  are  low  at  the  satellite  selection  time,  that 
become  very  large  after  the  selection  time.  Thus,  the  sat¬ 
ellite  set  has  larger  position  errors  over  the  interval  be¬ 
tween  satellite  selection  times.  This  effect  could  also 
occur  for  a  satellite  set  chosen  by  GDOP,  although  GDOP  has 
no  knowledge  of  the  measurement  noise  to  use  in  satellite 
selection. 

The  cost  criterion  also  was  significantly  affected  by 
the  number  of  satellites  visible  at  the  satellite  selection 
times.  In  all  five  missions  the  cost  criterion  did  not 


choose  an  alternate  satellite  set  (with  respect  to  the  GDOP 
selected  satellite  set)  when  five  satellites  were  visible, 
but  did  choose  an  alternate  set  when  six  satellites  were 
visible,  even  under  widely  varying  noise  conditions.  This 
result  can  be  explained  by  the  fact  that  5  feasible  satel¬ 
lites  give  the  cost  criterion  only  5  satellite  sets  to  choose 
from.  For  most  realistic  geometry  and  jamming  scenarios  for 
5  visible  satellites,  the  geometrical  affects  will  outweigh 
the  measurement  noise  and  filter  covariance  weights  of  the 
cost  criterion.  Thus,  for  5  visible  satellites  the  cost  cri¬ 
terion  will  tend  to  choose  the  same  satellite  set  as  the  GDOP 
criterion.  The  cost  criterion,  therefore,  should  be  used 
when  six  or  more  satellites  are  visible.  This  satellite 
visibility  problem  poses  a  constraint  on  when  missions  could 
be  flown  (in  order  to  achieve  the  smallest  possible  naviga¬ 
tion  errors)  unless  the  number  of  satellites  in  the  GPS  con¬ 
stellation  could  be  increased. 

Recommendations 

The  conclusions  of  this  report  indicate  that  short¬ 
comings  of  the  cost  criterion  (the  limitations  of  satellite 
visibility  and  the  rate  of  change  of  the  noise  density)  must 
be  addressed  before  it  could  significantly  impact  the  GDOP 
criterion  used  in  the  present  GPS  design.  Specifically, 
further  research  into  the  performance  of  the  cost  criterion 
in  a  rapidly  varying  noise  density  environment  is  needed. 

Such  research  would  greatly  aid  in  modifying  the  cost 


criterion  to  overcome  the  noise  density  problem,  or  aid  in 
development  of  another  satellite  selection  algorithm. 

This  study  has  developed  a  valuable  design  and  analysis 
tool  in  the  modified  Integrated  GPS/INS  Computer  Simulation. 
It  is  strongly  recommended  that  the  modified  IGI  be  used  to 
understand  the  satellite  selection  problem  thoroughly,  to 
aid  in  the  design  of  new  satellite  selection  algorithms, 
and  to  test  these  new  algorithms  against  realistic  aircraft 
dynamics  and  jamming  threats.  Requests  for  the  modified  IGI 
Computer  Simulator  code  should  be  addressed  to  the  Reference 
Systems  Branch  of  the  Avionics  Laboratory,  ( AFWAL/AAAN ) , 
Wright-Patterson  AFB,  Ohio  45433. 


Introduction 


The  following  sections  of  this  appendix  will  present 
the  satellite  visibility  study  completed  in  support  of  this 
report.  This  satellite  visibility  study  is  based  on  earlier 
satellite  visibility  studies  of  the  24  satellite  constella¬ 
tion  (7:28,  31-33)  and  the  18  satellite  constellation  (non- 
uniformly  spaced  satellites  within  each  of  the  3  orbits) 

(7:29,  34-36).  This  satellite  visibility  study  will  deter¬ 
mine  the  number  of  GPS  satellites  that  are  visible  to  the 
user  for  the  present  18  satellite  constellation  (uniformly 
spaced  satellites  within  each  of  the  6  orbits).  Since  the 
number  of  satellites  available  to  the  user  is  a  function  of 
the  user's  position  and  satellite  mask  angle,  these  parameters 
will  be  emphasized  in  the  satellite  visibility  discussion. 

Purpose 

The  purpose  for  a  satellite  visibility  analysis  as  part 
of  this  report  is  to  determine  if  the  number  of  satellites 
available  to  the  user  would  support  the  added  complexity  and 
computational  loading  of  an  alternate  satellite  selection 
criterion  (with  respect  to  the  standard  GDOP  satellite  selec¬ 
tion  criterion) .  The  number  of  visible  satellites  determines 
the  number  of  four  satellite  sets  from  which  the  best  sat¬ 
ellite  set  (with  respect  to  the  satellite  selection  criterion 
employed)  is  chosen.  If  there  is  a  large  number  of  satel¬ 
lites  visible,  there  is  a  very  large  number  of  satellite 


sets  available  and  the  large  number  of  computations  associ¬ 
ated  with  an  alternate  satellite  selection  criterion  may  be 
unacceptable  for  analysis  and/or  implementation.  If  there 
are  only  four  or  five  satellites  visible  the  small  number  of 
four  satellite  sets  (0  or  5  respectively)  effectively  pre¬ 
cludes  any  choice  of  an  alternate  satellite  set  regardless 
of  the  satellite  selection  criterion.  Also,  most  jammers 
will  be  either  ground-based  or  at  relatively  low  elevation 
angles  with  respect  to  the  GPS  satellites.  Therefore,  by 
masking  out  all  signals  at  low  elevation  angles  satellites 
with  high  noise  levels  (due  to  jamming)  and/or  deceptive 
jammers  (i.e.,  that  simulate  GPS  signals)  are  excluded  from 
satellite  selection.  Therefore,  satellite  visibility  is  of 
major  importance  in  the  selection  and  analysis  of  any  alter¬ 
nate  satellite  selection  criterion. 

Method 

The  method  used  to  determine  satellite  visibility  was 
a  computer  program  used  in  the  satellite  visibility  studies 
of  Reference  7.  The  computer  program  was  obtained  through 
the  Avionics  Laboratory's  Reference  Systems  Branch  (Air 
Force  Wright  Aeronautical  Laboratories,  Wright-Patterson  AFB, 
Ohio).  The  program  is  composed  of  software  modules  extracted 
from  the  Integrated  GPS/Inertial  Simulator  Computer  Program 
( IGI )  (17). 

The  major  modification  to  the  satellite  visibility 
computer  program  was  the  conversion  to  parameter  values  that 
describe  the  present  18  satellite  constellation.  The 


parameters  that  describe  the  18  satellite  constellation 


are: 

(1)  6  orbital  planes;  3  satellites  evenly  spaced 
within  each  plane 

(2)  circular  orbits  with  an  inclination  of  55° 

(3)  the  period  for  one  orbit  is  approximately  12  hours 

(4)  the  altitude  of  the  satellites  is  10,500  nautical 

miles 

(5)  the  orbital  phasing  is  40°  (in  an  eastward  direc¬ 
tion)  . 

In  addition  to  the  satellite  constellation  parameters 
the  desired  mask  angles  (there  are  8  possible  mask  angles 
that  can  be  specified  each  run)  were  varied  in  the  computer 
program. 

£ 

Satellite  Visibility  Analysis 

As  stated  in  the  introduction,  the  satellite  visibility 
is  dependent  on  the  user's  position  and  the  desired  mask 
angles.  The  satellite  visibility  discussion  will  look  at  the 
average  number  of  satellites  for  a  24  hour  observation  period 
for  latitudes  of  0°  to  90°  and  longitudes  of  10°  and  -90°. 

Case  1 .  The  satellite  visibility  for  a  user's  posi¬ 
tion  of  -90°  longitude  and  a  range  from  0°  to  90°  latitude 
(for  1°  increments  of  latitude)  is  examined  in  Figure  A-l. 
From  Figure  A-l  the  number  of  satellites  visible  for  mask 
angles  of  0°,  10°,  and  20°  is  approximately  7,  5.5,  and  4 
satellites  respectively  (for  0°  to  90°  latitude) .  The  sat- 
>V-  ellite  visibility  for  mask  angles  of  0°  to  20°  is  of 
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Figure  A-l .  Satellite  Visibility  -  55  Degree  Orbit 


particular  interest  in  determining  the  number  of  four  satel¬ 
lite  sets  available  to  a  satellite  selection  criterion. 

Figure  A-2  examines  satellite  visibility  for  mask  angles  be¬ 
tween  0°  and  25°.  Note  that  for  mask  angles  between  5°  and 
15°  the  average  number  of  visible  satellites  for  all  lati¬ 
tudes  is  between  5  and  6  satellites.  This  indicates  that 
mask  angles  for  10°  to  20°  are  undesirable  for  use  in  an 
alternate  satellite  selection  criterion  because  the  number 
of  satellite  sets  prevents  any  decisions  (other  than  geome¬ 
try  factors)  from  being  made  (See  the  Purpose  Section  of 
this  Appendix) .  Also,  for  mask  angles  between  5°  to  10° 
there  is  very  little  difference  in  the  number  of  visible  sat¬ 
ellites  for  increments  of  the  mask  angle.  Therefore,  since 
the  present  GPS  UE  uses  a  mask  angle  of  5°,  it  is  appropriate 
to  use  this  mask  angle  for  analysis  of  alternate  satellite 
selection  criterion. 

From  the  data  in  Figures  A-l  and  A-2  an  analysis  of  sat¬ 
ellite  visibility  for  a  user's  position  of  38°  Latitude, 

-90°  Longitude  is  possible.  This  location  is  of  particular 
importance  to  this  report  because  it  is  the  initial  user's 
position  in  the  F4  Close-Air-Support  flight  trajectory  used 
in  the  modified  IGI  Simulator  as  discussed  in  Chapters  IV, 

V,  and  Appendix  B.  The  number  of  visible  satellites  for  mask 
angles  from  0°  to  75°  is  shown  in  Table  A-l. 


18-Satellite  Uniform  Constellation 


Figure 


TABLE  A-l 


Happing  of  Satellite  Visibility  for  a  24  hr.  Period 
Latitude  38°  Longitude  -90° 


Mask  Angle  I  0°  5°  8°  10°  12°  15°  20°  25°  30°  45°  55°  65°  75° 


Av.  #  of  Sat. I  6.6  15.9  15.4  15.0  4.6  4.0 


3.4  3.0 


0.50.0 


Again  the  number  of  visible  satellites  for  mask  angles  from  10 
to  15°  would  result  in  the  average  number  of  satellite  sets  be 
tween  1  and  5,  undesirable  for  use  with  an  alternate  satel¬ 
lite  selection  criterion.  For  a  mask  angle  of  five  degrees 
the  user  could  anticipate  an  average  of  about  six  visible  sat¬ 
ellites  (15  possible  satellite  sets)  for  the  analysis  of  an 
alternate  satellite  selection  criterion.  With  15  satellite 
sets  the  alternate  satellite  selection  criterion  should  be 
effective  without  unacceptable  computational  loading. 

Case  2.  An  analysis  of  the  satellite  visibility  for  user 
positions  of  10°  Longitude  and  0°  to  90°  Latitude  is  displayed 
in  Figure  A-3 .  Figure  A-3  shows  some  variation  of  satellite 
visibility  from  that  shown  in  Figure  A-2,  but  the  number  of 
visible  satellites  is  consistent  with  the  values  seen  pre¬ 
viously.  A  tabular  view  of  the  number  of  visible  satellites 
for  mask  angles  from  0°  to  75°  is  shown  in  Table  A-2. 


TABLE  A-2 

Mapping  of  Satellite  Visibility  for  a  24  hr.  Period 
Latitudes  0°-90°  Longitude  10° 
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Comparisons  of  Table  A-l  and  A-2  show  a  slight  variation  in 
the  number  of  visible  satellites  with  varying  longitudes. 

The  variations  may  be  significant  for  a  particular  mission 
but  are  not  significant  for  the  analysis  of  alternate  satel¬ 
lite  selection  criterion. 

Summary 

The  results  of  a  satellite  visibility  study  have  been 
presented  in  this  Appendix.  The  purpose  for  this  study  was 
to  determine  if  an  alternate  satellite  selection  criterion 
would  be  acceptable  for  analysis  and/or  implementation  (in  an 
integrated  GPS/INS)  given  the  number  of  visible  satellites. 
Satellite  visibility  implies  two  constraints  on  an  alternate 
satellite  selection  criterion.  The  first  constraint  occurs 
when  there  are  four  or  five  satellites  visible.  This  small 
number  of  satellites  means  there  are  few  or  no  decisions  for 
the  alternate  satellite  selection  criterion  to  make  (other 
than  a  simple  geometry  decision) .  For  this  case  an  alternate 
satellite  selection  criterion  would  not  be  warranted.  The 
second  constraint  occurs  when  too  many  satellites  are  visible 
A  large  number  of  satellites  means  a  very  large  number  of 
possible  satellite  sets.  For  such  very  large  numbers  of  sat¬ 
ellite  sets,  the  increased  computational  loading  of  an  alter¬ 
nate  satellite  selection  criterion  would  be  unacceptable  for 
analysis  and  implementation. 

The  satellite  visibility  results  in  Case  1  and  2  show 


that  an  average  of  6  satellites  could  be  expected  for  a  mask 
angle  of  5°  (for  two  sets  of  user  position:  10°  Longitude 


0°  to  90°  Latitude  and  -90°  Longitude  0°  to  90°  Latitude). 

The  satellite  visibility  number  is  however,  an  average  over 
a  24  hour  observation  period,  this  means  that  periods  of 
less  than  6  satellite  visibility  occur,  and  that  during  these 
periods  the  alternate  satellite  selection  criterion  would 
probably  be  ineffective. 

Although  the  satellite  visibility  analysis  that  is 
presented  in  this  appendix  examines  only  two  sets  of  user 
positions  the  results  are  consistent  with  satellite  constel¬ 
lation  analyses  of  Reference  4. 


Appendix  B:  Modifications  to  the  IGI  Simulator 

The  following  sections  in  this  appendix  will  present 
the  assumptions,  algorithms,  and  equations  used  to  modify 
the  IGI  Computer  Program  Simulator.  The  modifications  fall 
into  three  main  areas;  the  null  steering  antenna,  adaptive 
bandwidth  control  of  the  satellite  tracking  loops,  and  the 
jamming  field.  The  implementation  of  the  cost  criterion 
(Equation  31)  will  not  be  discussed  in  this  appendix. 

Null  Steering  Antenna 

Introduction.  As  stated  in  Chapter  II’s  User  Segment 
discussion  the  null  steering  antenna  algorithm  acts  to  dis¬ 
tribute  its  nulls  to  minimize  the  total  power  out  of  the 
antenna  (6:139).  This  process  is  accomplished  by  combining 
the  outputs  of  several  individual  antenna  elements  through 
phase  and  amplitude  controls  to  minimize  the  power  out  of 
the  antenna  (10:6).  The  maximum  number  of  nulls  that  may  be 
formed  by  the  null  steering  antenna  is  dependent  on  the  num¬ 
ber  of  antenna  elements.  The  implemented  null  steering  an¬ 
tenna  model  uses  a  seven  element  phased  array  antenna  ar¬ 
ranged  in  a  circular  pattern  (Figure  B-l) .  This  seven  ele¬ 
ment,  circular  pattern  is  consistent  with  the  null  steering 
antenna  being  used  in  GPS  High  Performance  Aircraft  Develop¬ 
ment  and  Operational  Test  and  Evaluation  (25:1-3).  A  seven 
element  phased  array  antenna  may  create  a  maximum  of  six 
nulls  (6:139).  Each  of  the  antenna  elements  is  omnidirec¬ 
tional  in  azimuth  only.  The  antenna  elevation  gain  for  the 
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Figure  B-l .  Seven  Element  Null  Steering  Antenna  Geometry 


individual  antenna  elements  is  given  by  the  following 
equation  (6s 139) s 


for  0  <  80° 

g(0)  =4  cos  0 

(in  dB) 

(55a) 

for  0  >  80° 

g ( 0 )  =  1/3(90-0) 

-2.64  dB 

(55b) 

where  the  angle  0  represents  the  depression  from  the  normal 
(vertical)  to  the  aircraft.  Equation  55b)  accounts  for 
blockage  caused  by  the  aircraft  (6:139-140).  The  angular 
orientation  of  the  antenna  aperature  is  assumed  to  be  fixed 
with  respect  to  local  level  (i.e.,  the  aircraft  may  be 
oriented  other  than  straight  (the  pitch  axis)  and  level  (the 
roll  axis),  but  the  antenna  apperature  is  always  straight 
and  level  with  respect  to  the  local  level  coordinate  frame). 
This  assumption  was  made  to  simplify  the  software  to  describe 
the  null  steering  antennas  physical  characteristics.  A 
more  accurate  representation  of  the  physical  model  of  the 
antenna  will  require  coordinate  transformation  of  the  antenna 
elements'  position  from  aircraft  frame  to  local  level.  This 
transformation  must  be  calculated  using  filter  estimated 
angular  orientation.  The  null  steering  algorithm  would  then 
have  to  keep  track  of  measurement  noise  estimates  based  on 
the  true  antenna  elements'  position  and  the  filter  estimated 
antenna  elements'  position.  This  level  of  complexity  and 
impact  on  computer  processing  time  was  deemed  unappropriate 
for  this  analysis. 


Output  Pover  Minimization  Algorithm.  The  process  of 
minimizing  the  total  power  out  of  the  antenna  is  accomplished 
by  combining  the  outputs  of  several  individual  antenna  ele¬ 
ments  through  phase  and  amplitude  controls  (10:6).  This 
type  of  phase  and  amplitude  control  is  based  on  the  concept 
of  mutual  coupling  between  receiving  elements  in  a  phased 
array.  For  example,  two  elements  in  proximity  to  each  other 
induce  voltages  at  the  terminals  of  each  other  and  the  re¬ 
flected  voltage  in  each  line  depends  upon  the  phase  of  the 
driven  voltage  in  the  other.  Thus  in  a  phased  array  of  many 
elements,  the  reflection  coefficient  in  any  one  element  is 
unavoidably  influenced  by  mutual  coupling  and  the  total  ef¬ 
fect  will  depend  upon  the  phase  programming  of  all  the  ele¬ 
ments  in  the  phased  array  (26:1662).  The  phase  and  amplitude 
programming  is  accomplished  by  weighting  the  input  signal  to 
each  antenna  element.  The  input  signal  vector  (x)  to  the 
antenna  is  a  function  of  the  jamming  power  (J),  the  satellite 
signal  power  (S),  and  the  environment  noise  level  (a^)  at  the 
antenna  ( 6 : 140 ) . 


X  =  cr,  +  £  J  V  +  E  S.T  Y.t 
—  b  m  m  _m  n  N  N 


(56) 


The  vectors  and  are  complex  vectors  for  the  jammers 
and  satellites  respectively,  which  contain  the  phase  shift 
due  to  the  antenna  element  geometry  (6:140). 


V1  =  exp( il  • r 1 )  -  ifc  element  of  vector  V_ 
m  m  “ 


(57a) 


-v  '.  A V.'-', 


=  exp(j£n*r1)  =  1  element  of  vector  Y^ 


(57b) 


where: 

i  th 

r  =  position  of  the  i  omni  element  expressed  in  the 

local  level  coordinate  frame 

til 

1^  =  line-of -sight  to  the  m  jammer  expressed  in  the 

local  level  coordinate  frame 

til 

=  line-of -sight  to  the  n  satellite  expressed  in 
the  local  level  coordinate  frame 
The  background  noise  level,  cr^,  is  set  to  -l33dBw  and  the 
nominal  value  of  SN  =  -l62dBv  (6:140  and  5:148). 

The  input  signal  vector  X  (Equation  56)  may  now  be 
used  to  form  the  signal  autocorrelation  matrix  Rx  (6:104, 
27:81) . 


Rx  =  E(}QC*) 


(58a) 


R  =  a,  2 1  +  Z  V  *  +  E  S  2Y.tY.* 
— X  b  „  — m— m  KT  N  —N—N 

mm  N 


(58b) 


where  the  indicates  the  conjugate  transpose.  The  weights 

are  formed  from  the  inverse  of  Rx  and  a  control  law  de- 

A 

pendent  vector  e. 


-1  A 
W  =  Rx  e 


(59) 


£  =  [0,0,0,0,0,0,l]T 
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The  vector  e  imposes  the  constraint  that  one  of  the  elements 
have  unity  weight  and  the  remaining  elements  have  complex 
weights  (the  unity  weight  is  the  center  element  of  Figure 
13)  (6sl41).  For  computer  implementation  Equation  (59)  is 
rearranged  to  an  AX  =  B  form  for  solution  by  a  high  accur¬ 
acy,  complex  matrix,  linear  equation  solver  (28).  The 
weighting  vector  in  Equation  (59)  is  normalized  every  time 
the  null  steering  algorithm  is  called  to  prevent  degradation 
of  the  GPS  signals  due  to  old  weights  (serious  degradation 
of  the  GPS  signals  could  result  when  in  a  rapidly  varying 
jamming  environment  antenna  element  weights  are  used  for  a 
relatively  long  time  before  recalculation  (usually  to  reduce 
computational  loading) . 

Carrier-to-Noise  Power  Density  (CAlo)  Estimation 

The  Jammer  Power  Estimate.  The  jammer  power  estimate 
at  the  antenna  output  (J  )  is  a  function  of  the  jammer 
power  received  at  the  antenna  input  (JR)  and  the  gain  as¬ 
sociated  with  the  received  jammer  power  (G(0)). 


JR  + 

g  ( e ) 

(db) 

(61) 

£(J « 

+ 

)  (db) 

(62) 

m  m 

m 

|s<») 

PI** 

Vf-ll2  (db) 

(63) 

where: 


til 

J„  =  the  power  of  the  m  jammer  at  transmission 
m 

Ls  =  the  transmission  (space)  loss  of  the  radiated 
m 

jamming  signal. 

g(0)  =  as  defined  in  Equation  (55a  and  b) 

* 

W  =  the  conjugate  transpose  of  the  antenna  element 
weighting  vector 

V  =  the  complex  power  vector  for  the  m  jammers  (see 
Equation  (57a)). 

The  received  power  from  each  jammer,  JRm,  is  calculated 
using  the  following  equation  (20:6-7): 


JTmG(0,  6)A 
4ttR2 


(64) 


where : 

m  4-  V» 

JR  =  the  received  power  of  the  m  jammer  in  watts 

m  4-  Vj 

JT  =  the  transmitted  power  of  the  m  jammer  in  watts 
G(0, 6)  =  the  gain  of  the  transmitting  antenna  in  the 
direction  0,0  of  the  receiving  antenna  (unit 
less) 

A  =  the  effective  area  of  the  receiving  antenna  in 
square  feet 

R  =  the  transmitter  to  receiver  distance  in  feet. 
From  basic  antenna  theory,  the  effective  area  of  an  ideal 
isotropic  antenna  is  (20:7), 


A  =  X2/4tt 


(in  square  feet) 


(65) 


(66 


\  ^  =  wavelength  in  feet 

where: 

c  =  speed  of  an  electromagnetic  wave  in  a  vacuum  = 
9.8357119  x  108  feet 
f  =  frequency  of  transmitted  jammer  power,  assumed  to 
be  GPS  LI  =  1.57542  x  109  Hertz 
For  an  antenna  with  onmidirectional  gain  in  the  upper  hemi¬ 
sphere  the  gain  (G(0, 0)  is  a  constant.  The  constant  gain 
effect  permitted  the  use  of  Effective  Radiating  Power  (ERP) 
as  the  input  to  the  computer  simulation  of  the  jammer's 
transmitted  power  (20:7). 

ERP  =  vTTinG  < 0”,  © )  =  a  constant  input  by  the 

computer  simulation  user  (67 

Using  Equations  (65)  -  (67)  and  substituting  back  into  Equa 
tion  (64)  results  in  the  following: 

Jpm(dB)  =  ERP(dB)  +  [10  log( \2/(4n) 2 ) ]-20  log  R  (68a 

=  ERP(dB)  -  ( 26 . 076dB) -20  log  R  (68b 

Rewriting  Equation  (61)  by  substituting  in  Equations  (68b) 
and  (63)  results  in  the  following: 
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1  K  \9  »  .  r ,  ir^'r’’  r.  *  J  ^  A  •  ♦  .  •  .  \  ?  ■  ■  y  -tt ■  i  ^  ■  i  « i  1  i 


m 

£  (ERP. -(26 .076) -20  log  R  + 
i=l  1 


9^(9) 


(69) 


The  gain  associated  with  Equation  (63),  that  appears  as  the 
last  term  in  Equation  (69)  is  limited  to  a  maximum  null 
depth  of  30db  watts  (gain  =  30db  watts)  based  on  estimates 
of  broadband  jammer  power  (6:141). 

The  Signal  Power  Estimate.  The  signal  power  estimate 
at  the  antenna  output  (S  )  is  similar  to  the  approach  used 

cl 

in  the  jammer  power  estimate. 


Sa  =  SR  +  G<6> 

(dB) 

(70) 

SR  =  ST  +  LS  - 

2(9) 

(dB) 

(71) 

G ( 9 )  =  g( 9)  2  W 

*  2 

Y 

(dB) 

(72) 

where : 

ST  =  satellite  power  at  transmission 

=  ERP  (P-Code)  -  23.8db  watts  (Ref.  5:148) 

Ls  =  Space  Loss 

=  -101.75  -  20  log  R  (Ref  6:30) 

R  =  satellite  to  user  range 
2(9)  =  atmospheric  loss 

=  (0 .006 ) CC (6371  sin  9)2  +  127520]^-6371  sin  9]) 
(Ref  6:30) 

9  =  user  to  satellite  elevation  angle 
g(9)  =  as  defined  in  Equation  (55a  and  b) 
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Y^  =  the  complex  power  vector  for  GPS  satellite 
signal  (See  Equation  (57b)) 

The  space  loss,  Ls,  is  a  constant  -183.5db  for  most  near 
earth  missions  (Lg  varies  less  than  0.2  dB  for  user  altitudes 
up  to  200  miles).  The  atmospheric  loss  accounts  for  added 
satellite  attenuation  at  or  near  the  horizon.  This  term  can 
be  dropped  due  to  the  UE  set  masking  out  all  satellites  from 
the  horizon  to  5°  above  the  horizon  (i(0)  above  5°  from  the 
horizon  is  negligible  (less  than  -2db) ) .  Substituting  the 
above  simplifications  into  Equation  (71),  results  in  a  re¬ 
ceived  satellite  signal  of  -159. 7db,  slightly  better  than 
the  nominal  value  of  -162. Odb  used  in  Equation  (52b).  For 
simplicity  and  to  account  for  any  overly  optimistic  simplifica 
tions,  the  value  -162. Odb  was  used  for  SR .  The  resulting 
equation  for  the  signal  power  is. 


I 


Sa  =  (-162.0)  +  tg(0)l2LW*YNl2 


(73) 


'  *\ 


As  in  the  jammer  power  estimation  algorithm,  when  the  antenna 
places  a  null  on  top  of  a  GPS  signal  source,  the  null  depth 
is  limited  to  30db.  In  addition,  the  satellite  signal  gain 
is  increased  lOdb  in  the  range  -lOdb  <  G(0)  <  -30db  to 
simulate  the  effect  of  increased  gain  due  to  the  use  of 
directional  information  of  the  satellite's  position  in  the 
null  steering  controller  (reducing  the  chance  of  a  null 
being  put  on  a  satellite  signal,  thus  increasing  the  gain  on 
the  satellite  signal). 
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The  C/No  Estimate.  The  carrier-to-noise  density 
(C/No)  is  a  function  of  the  signal  power  estimate  (S  ),  the 
jamming  power  estimate  (J  ),  the  UE  set  front-end  noise 
power  density  (NT),  and  the  UE  set  processing  gain  (G)  as 
shown  in  the  following  equation  (5:150): 


P  s./K, 

=  - T - §— S - i -  (dB-H_) 

No  1  +  (Ja/sj(  a/M  )(1/r)  Z 


NT  7  v  ' G ' 


Nt  +  (Ja/G) 


(dB-Hz) 


where 

Nt  =  -202.0dbw/Hz  (Ref  5:148) 
G  =  70db-Hz  (P-code) 


(74a) 


(74b) 


The  term,  J  /G,  is  subject  to  the  assumptions  of  the  type  of 
jamming  signal.  In  this  implementation  it  is  assumed  that 
all  of  the  jamming  power  is  transmitted  into  the  satellite 
signal  20MHZ  bandwidth,  centered  at  the  frequency  of 
L2( 1.57542  x  109  HZ).  For  this  type  of  jamming  signal  the 
full  70db-HZ  of  the  P-code  processing  gain  is  used.  For 
jamming  signals  that  replicate  the  GPS  satellite  signal 
structure  the  processing  gain  may  be  reduced  significantly. 


Adaptive  Bandwidth  Control  of  the  Satellite  Tracking  Loops . 
Introduction .  The  purpose  of  incorporating  adaptive 


bandwidth  control  of  the  satellite  tracking  loops  is  to  simu 
late  the  contribution  this  algorithm  has  toward  the  GPS  UE 
Anti-Jam  capability,  simulate  the  two  modes  of  operation 


(coherent  and  non-coherent  tracking),  and  to  increase  the 
realism  and  capability  of  the  computer  simulation. 

The  Adapt ive  Bandwidth  Process .  In  the  User  Segment 
Section  the  importance  of  the  proper  satellite  tracking  loop 
bandwidth  on  anti -jam  performance  was  discussed.  The  adap¬ 
tive  bandwidth  process  attempts  to  minimize  the  mean  squared 
error  caused  by  all  disturbances  (noise,  jamming,  clock  in¬ 
stability,  and  aircraft  residual  dynamics)  affecting  the 
receiver  by  varying  bandwidth  as  a  function  of  C/No  (5:20). 

A  graphical  representation  of  the  adaptive  bandwidth  func¬ 
tions  for  the  coherent  and  non-coherent  mode  are  presented 
in  Figures  B-2a  dnd  B-2b  respectively.  For  implementation 
purposes  the  functions  were  simulated  by  two  second  order 
quadratic  equations  of  the  form 

Bc  =  A(C/Nq) 2  +  B(C/NQ)  +  C  (75) 

For  the  coherent  mode: 

For  ( 35 .Odb-HZ )  >  C/NQ  >  (29db-HZ) 

A  =  -0.0133 
B  =  1.1567 
C  =  21.20 

(max.  error  in  Bc  <  0.04HZ) 

For  (2 9. Odb-HZ)  >  C/Nq  >  (17db-HZ) 

A  =  0.0094 
B  =  0.3411 
C  =  3.1494 

(max.  error  in  Bc  <  0.07HZ) 


CODE  LOOP  NOISE  BANDWIDTH  (Hz) 


CARRIER  LOOP 


DELAY-LOCKED  LOOP 


ROUU* 


OH0  (dB-Hz) 


C/N0  (dB-Hz) 


Figure  B-2a  Bandwidth  Variation  with  C/JJ  for 
the  Coherent  Mode  (5:20)  ° 


SWITCH  MOOES 


SECOND  ORDER  LOOP 


DESIGN  POINT 
0.03  Hz  SdB  Hz 


C/Nq  (dB-Hz) 


Figure B-2b  Bandwidth  Variation  with  C/N  for 
the  noncoherent  Mode (5: 21)  ° 
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For  the  non-coherent  modes 

For  (17 .Odb-HZ )  >  C/NQ  >  (13db-HZ) 

A  =  0.0036  2  5 
B  =  0.0845 
C  =  0.563875 

(max.  error  in  B  <  0.0056HZ) 

c  — 

For  (13. Odb-HZ)  >  C/Nq  >  (6. Odb-HZ) 

A  =  0.00079167 
B  =  -0.007583 
C  =  0.04 

(max.  error  in  B  <  0.0021HZ) 

c  — 

From  the  figures  and  above  data  the  tracking  loop  limits  are 
readily  obvious .  Above  35dB-HZ  the  tracking  loops  are  oper¬ 
ating  in  an  environment  with  little  or  no  jammer  interference 
From  17.0  to  35.0dB-HZ  (and  above)  the  UE  is  in  the  coherent 
mode.  From  6.0  to  17.0dB-HZ  the  UE  set  is  in  the  non-coher¬ 
ent  mode.  Below  6.0dB-HZ  no  tracking  is  possible  (for  both 
satellites  being  tracked  or  for  acquisition  of  alternate  sat¬ 
ellites).  In  reality,  studies  have  shown  the  possibility  of 
tracking  a  satellite  below  6.0dB-HZ  until  correlation  errors 
grow  beyond  the  tracking  limits  (5:34-43).  The  modifications 
to  this  computer  simulation,  however,  implement  a  sharp  cut¬ 
off  at  each  lock  limit  to  deny  usage  of  satellites  with  mar¬ 
ginal  signal  strength  in  this  performance  analysis. 

Measurement  Error  Estimate .  The  measurement  error  for 
each  individual  satellite  measurement  is  directly  related  to 
the  RMS  error  of  the  tracking  loop  output.  The  RMS  error  of 


a  tracking  loop  output  is  composed  of  two  components:  0^ 
due  to  noise  and  jamming  (jitter),  and  o D  due  to  residual 
dynamics  and  oscillator  instabilities  (5:25).  The  total  RMS 
error  is  given  by: 

ae  *  (0N2  +  aD2)h  (76> 

The  jitter  term,  0^,  can  be  approximated  by  an  analytical 
expression.  The  0D  term  is  known  in  the  simulation,  but  is 
approximated  by  a  constant  1-0  value  in  the  IGI  Simulator. 

The  RMS  jitter  due  to  noise  and  jamming  for  the  coherent  mode 
can  be  approximated  by  the  following  expression  (5:25): 

0N  a  (98.424)  [l  +  37^]  %  (feet)  (77) 

where  B.  is  the  bandwidth  of  the  low  pass  filters  (these 
filters  separate  the  coded  signal  from  the  carrier  signal). 
The  expression  for  the  non-coherent  mode  is  given  by  multi¬ 
plying  the  (Bc/(C/toQ) )  term  in  Equation  (77)  by  0.5.  For 
computer  implementation  the  measurement  error  estimate  (Rm) 
is  approximated  by  the  RMS  error  of  the  track  loop  (0e) . 

The  one  sigma  truth  and  filter  measurement  errors  used  in 
the  IGI  simulator  are  7.0  feet  and  12.0  feet  respectively 
(for  an  environment  without  jamming) .  These  one  sigma 
errors  are  associated  with  the  0D  term  in  Equation  (76). 

The  one  sigma  filter  measurement  error  is  larger  than  the 
true  error  to  account  for  measurement  errors  associated  with 
atmospheric  disturbances  and  inaccuracies  in  the  satellites 


emphemeries  that  are  represented  in  the  truth  model  as  one 
sigma  position  errors  of  10  feet  (truth  states  58-61),  but 
are  not  modeled  in  the  filter  state.  Therefore  the  equa¬ 
tion  for  the  filter  measurement  error  (R  .)  is, 

mr 

Rjnf  =  (aN2  +  cD2  +  95.0)^  (feet)  (78) 

As  stated  in  the  previous  section  the  code  tracking  loop 
lock  limit  is  a  C/Nq  value  equal  to  6.0dB-HZ.  The  measure¬ 
ment  error  estimate  for  any  satellite  that  is  below  this 

Q 

lock  limit  is  set  to  a  value  of  1.0  x  10  feet.  This  measure 
ment  error  value  is  at  least  five  orders  of  magnitude  higher 
than  the  largest  filter  measurement  error,  which  should 
adequately  simulate  signal  loss-of-lock  and  appropriately 
deweight  the  satellite's  pseudo-range  and  delta  pseudo-range 
measurements  during  a  filter  update. 

The  Jamming  Field 

Introduction.  The  IGI  simulator  does  not  model  any 
jamming  field  related  parameters.  Therefore,  the  modifica¬ 
tions  to  the  simulator  could  exercise  a  great  deal  of  flexi¬ 
bility  in  modeling  the  jamming  field.  The  goal  behind  the 
modifications  to  be  discussed  in  the  next  subsection  is  to 
add  a  capability  that  is  consistent  with  this  analysis,  but 
capable  of  supporting  future  analysis  (using  realistic 
jamming  scenarios)  without  major  modifications. 

Jamming  Field  Parameters .  The  jamming  field  parameters 
added  to  the  IGI  Simulator  are: 


(1)  up  to  25  jammers  can  be  transmitting  at  one  time 

(2)  the  jammers  may  be  on  the  ground  or  airborne 

(3)  the  jammers  may  be  moving  or  stationary. 

Each  jammer’s  position  is  specified  in  degrees  of  latitude 
and  longitude,  and  altitude  in  feet.  Each  jammer’s  velocity 
is  specified  in  units  of  f eet-per-second  in  the  three  com¬ 
ponents  of  the  local  level  coordinate  frame  (i.e.,  Vjeast» 
Vjnorth*  and  V jup^  *  Power  (EIRP  or  sometimes  just  ERP)  is 
specified  in  kilowatts.  The  jammer's  have  an  omni-directional 
radiating  pattern  in  both  azimuth  and  elevation.  The  geo¬ 
graphical  terrain  (i.e.,  mountains,  valleys,  foliage)  that 
might  obstruct  the  user- jammer  line-of -sight  (LOS)  is  not 
modeled.  Instead,  the  LOS  limit  from  the  user-to- jammer  is 
determined  by  the  following  equation  (20:9): 

LOS  LIMIT  =  (8.0  R/3.0)^(  (H^  +  (H2)^)  (79) 

where : 

R  =  Radius  of  the  earth  (2.09  x  107  feet) 

H1  =  Jammer  height  (in  feet 

H2  =  The  user's  true  height  (in  feet) 

All  of  the  specified  jammer  parameters  are  listed  in 
the  initialization  portion  of  the  IGI's  printed  output. 

This  appendix  has  presented  the  assumptions,  algorithms, 
and  equations  associated  with  the  three  major  modifications 
to  the  IGI  Simulator.  Of  the  three  modifications  the  null 
steering  antenna  algorithm  was  the  most  difficult  to  imple¬ 
ment  in  software. 
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